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ABSTRACT 

Multifrequency observations from the VLA, VLBA and OVRO Millimeter Array of a major radio 
outburst of Cygnus X-3 in 2001 September are presented, measuring the evolution of the spectrum of 
the source over three decades in frequency, over a period of six days. Following the peak of the flare, 
as the intensity declines the high-frequency spectrum at frequency v steepens from v~ aA to ^ _0 6 , after 
which the spectral index remains at this latter terminal value; a trend previously observed but hitherto 
not satisfactorily explained. VLBA observations, for the first time, track over several days the expansion 
of a sequence of knots whose initial diameters are ~ 8 milliarcseconds. The light-crossing time within 
these plasmons is of the same order as the time-scale over which the spectrum is observed to evolve. 
We contend that properly accounting for light-travel time effects in and between plasmons which are 
initially optically thick, but which after expansion become optically thin, explains the key features of 
the spectral evolution, for example the observed timescale. Using the VLBA images, we have directly 
measured for the first time the proper motions of individual knots, analysis of which shows a two-sided 
jet whose axis is precessing. The best-fit jet speed is f3 ~ 0.63 and the precession period is ~ 5 days, 
significantly lower than fitted for a previous flare. Extrapolation of the positions of the knots measured 
by the VLBA back to zero-separation shows this to occur approximately 2.5 days after the detection of 
the rise in flux density of Cygnus X-3. 

Subject headings: radio continuum: stars — techniques: high angular resolution — radiation 

mechanisms: non-thermal — binaries: close — stars: individual (Cygnus X-3) 



1. INTRODUCTION 

Cygnus X-3 was first detected by Giacconi et al. (1967) 
using X-ray proportional counters during a 1966 rocket 
flight and rose to prominence with the detection of a giant 
radio flare in 1972 September (Gregory et al. 1972a). It 
is a high-mass X-ray binary system located in the Galac- 
tic plane at a distance of wlO kpc from the Earth (Pre- 
dehl et al. 2000; Dickey 1983). With Galactic co-ordinates 
I = 79.84°, b = 0.69°, it lies in or behind either the Perseus 
or the Outer Arm, depending on its exact distance from 
us. The nature of the compact object is still uncertain. 
The large interstellar extinction to the source precludes 
optical spectroscopy, making it difficult to obtain a re- 
liable mass function and hindering identification of the 
companion star. However, infrared spectra show strong, 
broad helium emission lines, and a lack of strong hydro- 
gen emission. Spectral analysis by van Kerkwijk et al. 
(1996) indicated that the companion is a Wolf-Rayet star 
of the WN7 subclass, although more recent observations 
by Koch-Miramond et al. (2002) suggested that the sub- 
class may in fact be WN8. The orbital period, as inferred 



from X-ray and infrared flux modulations, is 4.8 hours 
(e.g. Parsignault et al. 1972). This source occasionally 
undergoes huge radio outbursts where the flux density can 
increase up to levels of ~ 20 Jy. Of order 50 outbursts with 
peak radio flux densities exceeding 1 Jy have been observed 
since the detection of the first such flare in 1972 Septem- 
ber. During one of these radio outbursts, the hard X-ray 
emission was found to correlate strongly with the radio 
emission, although the soft X-ray emission was anticorre- 
lated with radio flux density (McCollough et al. 1999). In 
quiescence, the soft-X-ray is found to be strongly corre- 
lated with the radio emission, and anticorrelated with the 
hard X-rays (Choudhury et al. 2002). This soft-X-ray— 
radio correlation disappears immediately before and dur- 
ing periods of radio flaring (Watanabe et al. 1994). During 
recent outbursts, jet-like structures have been observed at 
radio frequencies, although there has been some debate as 
to the morphology of the collimated emission. Two-sided 
jets from a flare in 2000 September were detected on arc- 
second scales in a north-south orientation with the Very 
Large Array (VLA) (Marti et al. 2001), with an inferred 
jet speed v — 0.48c. On milliarcsecond scales, a one-sided 
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jet with the same orientation was imaged with the Very 
Long Baseline Array (VLBA) after a flare in 1997 Febru- 
ary (Mioduszewski et al. 2001). In this case, the derived 
jet speed was significantly higher (v > 0.81c). 

In this paper we analyse radio observations of the 2001 
September outburst of Cygnus X-3 and describe some of 
the physical mechanisms responsible for the observed evo- 
lution of the source. In §2 we summarize the characteris- 
tics of previous radio flares of Cygnus X-3. In §3, §4 and 
§5, we present our observations from the VLA, Owens Val- 
ley Radio Observatory (OVRO) Millimeter Array, and the 
VLBA respectively, and in §6 we analyse the VLBA images 
and infer the speed, sidedness and precession parameters 
of the jet. We go on to discuss the spectral evolution and 
the underlying physical mechanisms for the low-frequency 
turnover in §7. We then derive constraints on the magnetic 
field strength, electron number density and energy density 
in the source in §8, and describe the physics behind the 
high-frequency spectral evolution in §9. In §10 the mag- 
netic field strength and energy in the jet in this outburst of 
Cygnus X-3 are compared to derived parameters of radio 
galaxies such as Cygnus A and M 87. 

2. SPECTRAL CHARACTERISTICS 

The temporal evolution of the radio spectrum of Cygnus 
X-3 was first studied during the outburst of 1972 Septem- 
ber (Gregory et al. 1972a, and 20 papers following this). 
This outburst was observed between 0.4 GHz and 90 GHz, 
and observations were made during both the increase and 
the subsequent decrease of the radio flux density. This 
outburst of Cygnus X-3 exhibited approximately similar 
characteristics to several other major flares observed from 
this system: a short-lived exponential fading of the inten- 
sity (e.g. Aller & Dent 1972; Hjellming et al. 1974; Marsh 
et al. 1974; Waltman et al. 1995) with e-folding times r in 
the range 0.15 < r < 2.75 days, followed by a power-law 
decay in intensity. With sufficient frequency coverage a 
two-phase spectral evolution is seen; an initial steepening 
of the spectrum, after which the spectral index a (which 
relates the intensity S„ to the frequency v as S u oc v~ a ) 
remains at a terminal value of ~ 0.6 — 0.7 (e.g. Gregory 
et al. 1972b; Aller & Dent 1972; Hjellming & Balick 1972; 
Seaquist et al. 1974; Gcldzahler et al. 1983; Fender et al. 
1997). In the dataset we present in this paper, we have 
observations over a wider frequency range than made pre- 
viously, which more tightly constrain the observed evolu- 
tion. However, our observations were not triggered early 
enough to observe the rise phase of the outburst, which 
precludes detailed modelling of the mechanism responsi- 
ble for launching the jet, and comparison with previous 
models purporting to explain the rise phase of radio flares 
(e.g. Marti et al. 1992). 

3. VLA OBSERVATIONS AND DATA REDUCTION 

The outburst of Cygnus X-3 was observed to begin when 
the radio flux density started to rise on 2001 September 14 
(MJD 52167 ± 0.2) (Trushkin et al. 2002), and was mon- 
itored regularly with the VLA from September 18 until 
ten days following the outburst, and then less frequently 
until the end of October. The VLA was in its most com- 
pact "D" configuration, so that the source was unresolved 
and, crucially, the VLA measured the integral of all emis- 



sion regions associated with the outburst. The source was 
observed using snapshot observations of duration two to 
three minutes, sufficient to detect and perform photometry 
on the source. 

Data were taken with the VLA using standard proce- 
dures, observing at eight different frequency bands (cen- 
tred on 73.8 MHz, 327.5 MHz, 1.425 GHz, 4.86 GHz, 8.46 GHz, 
14.94 GHz, 22.46 GHz, and 43.34 GHz). A summary of the 
observing dates and frequencies may be found in Table 1. 
All observations except for those at 73.8-MHz were made 
in two independent frequency bands (IF pairs). All IFs 
were of width 50 MHz except for the 327.5- and 73.8-MHz 
data, which had bandwidths of 3.125 and 1.563 MHz re- 
spectively. At 1.425 GHz, the two IFs (with frequencies 
of 1.465 and 1.385 GHz for the upper and lower sidebands 
respectively) were imaged separately, to gain an extra low- 
frequency point where the spectral shape was changing 
most rapidly. At all other frequencies, the data in the sep- 
arate IFs were averaged together during imaging. Calibra- 
tion and image processing were performed using NRAO's 
Astronomical Image Processing System (AIPS). The pri- 
mary calibrators used were 3C 48 and 3C 286, depending 
on which was above the horizon at the time of the obser- 
vation. One of these was observed at each frequency in 
every observing run. The flux scale used was that derived 
at the VLA in 1999, as implemented in the 31Dec02 ver- 
sion of AIPS. Several different secondary calibrators were 
used, depending on the frequency of observation and on 
the other sources observed during the run. These were 
J 2007+4029, J 2015+3710, and J 1953+3537, at angular 
separations of 4.7°, 5.0°, and 9.3° from the source respec- 
tively. After initial calibration, the data were imaged and 
self-calibrated (initially using phase-only self-calibration, 
but ultimately using simultaneous amplitude and phase 
self-calibration to make the best possible images) in order 
that photometry could be performed on the source. 

3.1. Low frequency issues 

At low frequencies, the large antenna primary beam re- 
quired a slightly different observing strategy. The power- 
ful radio galaxy Cygnus A is located only 6.3° away from 
Cygnus X-3 and, at 17kJy, is the brightest source in the 
74-MHz sky. At 330 MHz, the flux density of Cygnus A 
is 5.3 kJy Different methods were employed at the two 
different frequencies. At 74 MHz, where the field of view 
is largest, a single observation was taken, with the point- 
ing centre located halfway between Cygnus A and Cygnus 
X-3. At 330 MHz, alternate scans were made on Cygnus 
A and on Cygnus X-3 as for the higher-frequency obser- 
vations. It was thought that the former method might 
allow the deconvolution process to remove the sidelobes 
from Cygnus A more accurately. The results suggested 
the latter method to be preferable, since a source was de- 
tected at the known position of Cygnus X-3 at 330 MHz. 
At 74 MHz however, it was only possible to put an upper 
limit on the flux density of the source, taken as three times 
the r.m.s. background noise level of the final sky image at 
the location of Cygnus X-3. 

3.2. High-frequency issues 
3.2.1. Opacity corrections 
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At the highest frequencies (43 GHz and in some cases, 
22 GHz), the weather had a significant effect on the data, 
causing poor phase correlation. This required scalar aver- 
aging or vector averaging over a short timescale (of order 
the integration time of 3.3 s) of the visibilities during the 
calibration process, rather than scan-based vector averag- 
ing as is customary. In addition, the variable atmospheric 
opacity required the use of weather (WX) tables or tip- 
ping scans (in which the system temperature for an indi- 
vidual antenna is measured as the antenna slews from the 
zenith to the horizon; the dependence of system tempera- 
ture on zenith angle may be used to determine the opac- 
ity) in order to correct for the effects of the opacity, and 
thus determine the flux scale accurately. Tipping scans 
were only available for the 22-GHz observations made on 
September 18 and 19. For all other epochs and frequen- 
cies, an estimate of the opacity had to be made using the 
weather information taken during the observations. The 
fitted opacity curves given by Butler (2001) were initially 
used for this. However, more recent analysis by Butler 
(2002) showed that an improved estimate of the atmo- 
spheric opacity may be made by including a term to ac- 
count for seasonal variations, which we therefore factored 
into the calculations. The use of a constant estimated 
opacity throughout an individual observation was justified 
by the short lengths of the observations. 

3.2.2. Referenced pointing 

Another problem in making high-frequency observations 
is the pointing accuracy of the antennas: the difference 
between the centre of the primary beam and the desired 
source position. Factors which may contribute to the point- 
ing error include atmospheric refraction, gravitational de- 
formation and differential heating of the antennas and 
non-perpendicularity or misalignment of the antenna axes 
(Clark 1973). Under normal observing conditions, the 
VLA pointing is accurate to 10-20", but can be as bad 
as an arcminute in poor weather conditions (Rupen 1997). 
At high frequencies, the primary beam is small (the full- 
width at half power being about 1' at 43 GHz and 2' at 
22 GHz). Thus the pointing error can be a large fraction 
of the primary beam, leading to an error in the observed 
amplitude of the source. To account for this, primary ref- 
erenced pointing was used in the observations (deriving 
the pointing offsets by pointing up on the calibrator at 
a lower frequency — typically 8.4 GHz — and applying 
the derived corrections). This lowers the r.m.s. pointing 
error to about 2" in azimuth and 5" in elevation at ele- 
vations below about 80°. Wind loading of the antennas 
may also affect the pointing, but the effect of winds be- 
low 6 ms -1 has been measured to be insignificant (Rupen 
1997). Fortunately, the maximum wind speed in any of 
our observations was 5.7 ms -1 . For the observing run early 
on September 21, the pointing accuracy in azimuth for a 
subset of the VLA antennas was poor at elevations > 80°, 
leading to significant gain fluctuations. After removing the 
affected data, we found that the point-source flux densi- 
ties before and after amplitude self-calibration agreed to 
within 1%. 

3.2.3. Gain corrections 

As the VLA antennas track in elevation, the variation 
in gravitational deformation of the antennas causes the 



reflector profiles to change. This makes a difference to the 
response of the telescope at high frequencies (at and above 
15 GHz, where the wavelength is more comparable to the 
scale of the distortions). To account for this measurable 
effect, the gain curves closest in time to the observing dates 
(those of 2001 November), as measured by NRAO staff, 
were used. 

3.3. Errors and Uncertainties 

There are several possible sources of uncertainty in the 
data reduction process which must be taken into account 
when estimating the uncertainties on the VLA measure- 
ments. The r.m.s. noise in the final image of the source 
gives an idea of the random error in the measurements 
due to thermal noise fluctuations. For 2-minute integra- 
tions with the full VLA, the expected thermal noise is at 
best 0.1 mJy per beam at 8.4 GHz. At higher frequen- 
cies, the system temperature increases due to atmospheric 
emission, so the noise rises to 0.6 mJy per beam at 43 GHz. 
At low frequencies, the Galactic and cosmic backgrounds 
increase the system temperature markedly, to 34 mJy per 
beam at 330 MHz and 3.3 Jy per beam at 74 MHz for ob- 
servations in the Galactic plane. In practice, even this sen- 
sitivity level will not be reached for low-frequency snapshot 
imaging, since in this case it is impossible to solve accu- 
rately for the sidelobes of confusing sources with only a 
short observation. 

In addition to the thermal noise in the image, there are 
several sources of uncertainty resulting from the calibra- 
tion process, the most significant of which is likely to be 
the uncertainty in the determination of the flux scale. The 
AIPS tasks SETJY and GETJY quote nominal uncertain- 
ties in the calibrator fluxes, but these are based solely on 
the scatter in the amplitude gain ratios. If there are other 
effects, such as phase decorrelation of a weak phase calibra- 
tor (such as J 1953+3537 at high frequencies, which was 
used occasionally) or due to bad weather, then depending 
on the form of averaging and the solution interval used in 
the AIPS task CALIB, the derived antenna voltage gain 
solutions and hence the flux scale will vary. 

CALIB derives the complex antenna voltage gain so- 
lutions for amplitude and phase for sources with known 
flux densities and structures, thus setting the amplitude 
scale. To reduce the error in the derived gains, the cali- 
brator data are averaged in time. To obtain an unbiased 
estimate of the phase, vector averaging is used. However, 
it is more difficult to obtain an unbiased estimate of the 
amplitude, since vector averaging will potentially underes- 
timate the amplitude due to decorrelation by phase noise. 
Alternatively, scalar averaging will tend to increase the 
amplitude due to the inclusion of a noise bias; a scalar 
average of a data set consisting of pure noise will pro- 
duce a positive signal. However, this is only a serious 
problem for high-frequency observations, where the phase 
decorrelation tends to be greater, and for weak (< 5a) 
sources. During the observations in September, Cygnus 
X-3 was bright (w 10 2 — 10 3 cr), so scalar averaging should 
not introduce a significant noise bias. Since the amplitude 
and phase fluctuations of most instrumental effects are un- 
correlated (Fomalont & Perley 1999), separate amplitude 
and phase averaging may be used to give the most accu- 
rate estimates of each (although when deriving the actual 
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antenna-based solutions for the complex gains, the noise 
bias is reduced by solving simultaneously for both). In 
the data reduction process, the initial approach was to use 
vector averaging over a whole scan, solving for amplitude 
and phase simultaneously. At high frequencies (> 15 GHz) 
scalar averaging was then also done as a check. If the two 
disagreed, separate calibration was done for phase and am- 
plitude, solving for phase first using vector averaging, and 
then applying those gains and solving simultaneously for 
amplitude and phase using scalar averaging. 

According to the calibration procedure laid out in the 
VLA calibration manual, the flux density bootstrapping is 
accurate to 1-2% for the 1.4- , 4.8-, and 8.4-GHz observing 
bands. The uncertainty in the measurements was taken 
as 1% for 8.4 GHz and 2% for 1.4 and 4.8 GHz. At higher 
observing frequencies, the accuracy is possibly as good as 
3-5%, but effects such as poor antenna pointing and atmo- 
spheric opacity are thought to reduce the accuracy further. 
From the data reduction, it was clear that the atmosphere 
was particularly troublesome on certain days (September 
18, 21, 23, and 24), and phase decorrelation made the 
flux scale uncertain. On those days, the percentage uncer- 
tainty was taken as the percentage flux density difference 
between the GETJY flux densities of the phase calibrators 
found using vector averaging and scalar averaging during 
the calibration. This difference was only found to be sig- 
nificant at 43 GHz. For the better days (September 19 and 
22), and for the observations at 15 and 22 GHz, the uncer- 
tainty was derived by spectral fitting. These uncertainties 
were then combined in quadrature with the thermal noise 
in the images (although since we achieved a dynamic range 
of several hundred to one in most of the images, the con- 
tribution of the thermal noise to the total uncertainty was 
minimal). 

It seems reasonable to assume that the spectrum is smooth. 
Fitting a power law to the spectrum at frequencies well 
above the low-frequency turnover, the deviations of the 
high-frequency points from the expected smooth curve were 
measured. The uncertainty in the measurement for a par- 
ticular frequency was then taken as the r.m.s. deviation 
of the points from their respective power laws (taking the 
mean over all six days). This resulted in uncertainties of 
4.8% for the 15-GHz data, 7.7% for the 22-GHz data, and 
8.0% for the 43-GHz points not affected by phase decor- 
relation (for the data so affected, the uncertainty was of 
order 7-12%). 

At low frequencies (74 MHz and 330 MHz), the flux scale 
was determined from observations of Cygnus A. It is thought 
to be good to 5-10% (R.A. Perley 2003, private communi- 
cation) . 

4. OVRO OBSERVATIONS AND DATA REDUCTION 

Five observations of Cygnus X-3 were made by the OVRO 
Millimeter Array during the period September 20 to 26 
(see Table 2), each involving between 1 and 3 hours of in- 
tegration. This array is comprised of six 10.4 m antennas 
on a 400 m T-shaped track, situated at 1220 m altitude in 
the high desert of eastern California. The array was in 
compact configuration during these observations, and had 
a typical angular resolution of ~ 11 "x 7.5". Two 1-GHz 
sidebands centred at 97.98 and 100.98 GHz were observed 
at a single polarisation at all epochs. The flux density scale 
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Fig. 1. — VLA and OVRO Millimeter Array observations of the 
integrated flux density during the outburst of Cygnus X-3. The bold 
black horizontal lines indicate the times of the VLBA observations. 
The arrow indicates the start of the outburst, as detected with the 
RATAN-600 telescope (Trushkin et al. 2002). Lines have been drawn 
between the points corresponding to September 9 and September 18, 
to indicate the change in spectral shape. This does not represent the 
rate at which the flux density rose during the outburst, as there is 
insufficient temporal coverage to accurately plot the rise phase of the 
flare. Note how the spectrum goes from rising with frequency before 
the flare to falling with frequency afterwards, a trend previously 
observed as the source moves from the prc-flarc "quench" state to 
outburst (e.g. Fender et al. 1997). A power-law fit to the high 
frequency decay curves between September 19 and September 23 
yields a power-law index of 3.8 ± 0.2. 



was set by assuming the flux density of the adjacent con- 
tinuum calibrator J 2015+3710 was 2.87 Jy for all epochs 
(this value was derived from the first observations made 
on September 20 involving Uranus and 3C345, and incor- 
porates the on-line opacity correction), and is considered 
accurate at the 10% level. Reduction of the data was per- 
formed using the NRAO AIPS package. The typical la 
image detection limit for the observations was 5-10 mJy. 
The lightcurves obtained using the VLA and OVRO Mil- 
limeter Array are shown in Figure 1 . 

5. VLBA OBSERVATIONS AND DATA REDUCTION 

Cygnus X-3 was also monitored with the VLBA during 
the September outburst. Each day from 2001 September 
18 to 23, observations of duration between 8 and 13 hours 
were made using all ten VLBA antennas and a single dish 
of the VLA. The observation periods are shown in Figure 
1. On all six days, observations were made at 1.660, 4.995, 
15.365, and 22.233 GHz, with the time being split approxi- 
mately equally between the different frequencies. Scans at 
a given frequency were spread over the whole observation 
period in order to improve the ww-coverage. 

The VLBA observations were carried out using dual po- 
larisation and two-bit sampling. At all four frequencies, 
an 8-MHz bandwidth was used. The two IF pairs were 
allocated across an uninterrupted frequency range, rather 
than being more widely spaced. Scans on Cygnus X-3 were 
interleaved with scans on a phase calibrator source, with 
a 3-minute cycle time (70s on the calibrator, 110s on the 
source). The calibrators used were J 2052+3635 (5.9° from 
the source) for September 18 and 19, and J 2007+4029 
(4.7° from the source) for all other epochs. The data were 
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correlated using the VLBA correlator in Socorro, New 
Mexico. All data taken at <15° elevation were nagged. 
Amplitude calibration, fringe-fitting and imaging were car- 
ried out using the AIPS software according to NRAO 
guidelines laid out in Appendix C of the AIPS Cook- 
book. The flux density scale was set using system tem- 
peratures and antenna gains, and is believed to be accu- 
rate to <~ 5% according to the VLBA Observational Status 
Summary. Fringe-fitting was carried out only on the cali- 
brator sources, not the source itself, except at 22 GHz as 
described later in this section. 

Unfortunately, Cygnus X-3 is one of the most heavily 
scatter-broadened sources in the radio sky owing to its 
location in the direction of the Cygnus OB2 association 
(Wilkinson et al. 1994; Mioduszewski et al. 2001). The size 
of the scattering disk scales as v~ 2 , so particularly at the 
lower observing frequencies, much of the information from 
the long baselines was lost, and the resulting uu-coverage 
of the observations was poor. In addition, the images of 
the phase calibrators themselves were heavily scattered, 
so it was not possible to calibrate the outer antennas at 
low frequencies. Only the baselines between the anten- 
nas in the southwestern United States (the VLA antenna, 
Pie Town, Los Alamos and Fort Davies) yielded useful 
information about the source. This meant that we were 
unable to produce any good quality images at 1.6 GHz, 
and that the resolution decreased more markedly with fre- 
quency than would have been expected from a simple X/D 
decline. 

In addition to sparse wu-coverage (the VLBA itself has 
a maximum of ten antennas, the outermost of which were 
lost owing to the scattering of information from long base- 
lines), the complex, variable nature of the source during 
each observation posed further difficulties for the data re- 
duction. In an observation of duration 13 hours, the com- 
ponents of the source would not only change in flux den- 
sity, but would also gradually change position, owing to 
the motion of the jet during the period of observation. 
Modelling the latter effect however revealed that the po- 
sitional changes would simply elongate the components in 
the images. This effect would in fact only be a problem at 
the highest frequency, since for a 13-hour observation, a 
proper motion of ~ 10 milliarcseconds (mas) per day would 
shift the position by 6.5 mas. This is a little smaller than 
the beam size at 22 GHz (see Table 3), significantly less 
than the beam size at 15 GHz and much smaller than the 
size of the scattering disk at 5 GHz. So even at 22 GHz, the 
effect would be fairly small. The earlier epochs of VLBA 
data suffered most from source variability, whereas the rel- 
atively low surface brightness of the source caused most of 
the problems during the later observations. On September 
22 and 23, we were unable to produce any good images. 

Several approaches to producing good quality images 
of the source were tried. Splitting the observation into 
small timeranges and imaging each individually in order 
to reduce the effects of variability and proper motion was 
found to help for the early observations when the source 
was bright, but for the later epochs, did not leave enough 
flux density to make a good quality image of the source. 
Self-calibration with a point source model was attempted, 
but this was not found to improve the image significantly. 
Imaging only the part of the observation where a plot of 



the amplitude of the short baselines against time was rela- 
tively flat, in order to reduce the effects of variability, was 
again not found to improve the image. Downwcighting 
the central regions of the ira-plane (corresponding to the 
short baselines) by the use of super-uniform weighting was 
similarly found to be of limited use, as was the attempt 
to model the visibilities with Gaussian components. The 
method that did most to improve the map was to image 
and self-calibrate the phase calibrator (J 2007+4029), and 
then to transfer the solutions for amplitudes and phases to 
the target. Once the solutions had been transferred, some 
of the methods listed above were applied once more in 
an attempt to improve the images. Of these, only down- 
weighting the short baselines was found to significantly 
improve the images. The resulting images are shown in 
Figure 2. More information on the images (beam sizes, 
noise levels and exact observation times) is given in Table 
3. Only phase self-calibration was used since amplitude 
self-calibration was not found to improve the images. 

The observations at 5 and 15 GHz were phase-referenced 
to a calibrator (J 2052+3635 or J 2007+4029 as detailed 
above) in order to preserve positional information. Since 
self-calibration was carried out on both the phase calibra- 
tor and Cygnus X-3 itself, there might still have been some 
degree of uncertainty in the absolute sky positions of the 
images. In order to estimate this effect, a check source 
was imaged, both phase-referenced to the phase calibrator 
and to itself. The difference in the position of the check 
source in the two images gave an estimate of the positional 
uncertainty introduced in the imaging and self-calibration 
processes. The check sources used were J 2025+3343 for 
September 18 and 19 and J 2052+3635 for all other obser- 
vations. 

At 22 GHz, the phase calibrator sources were insuffi- 
ciently bright, and the images could not be phase- referenced 
Thus their absolute sky positions are not accurate; an in- 
evitable consequence of the fringing process is that the 
brightest point in the source becomes the centre of the 
image. At 15 GHz and 5 GHz however, phase referencing 
was carried out. 
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Fig. 2. — VLBA images of Cygnus X-3 in flare. Images are shown at 4.995 GHz in the left column, 15.365 GHz in the central column, and 
22.233 GHz in the right column. Where images are missing, the data were not good enough to make a good quality image. Time progresses 
down each column, each row a day later than the previous one. Solid contour levels are (\/2) n times the lowest level stated in the images 
(n = 0, 1, 2, . . .). Dashed contours are —1 and — v2 times the stated lowest levels. The beam sizes arc shown in the lower left corner of the 
images. Crosses indicate the fitted positions of the maxima, as plotted in Figure 4. Assuming a distance of 10 kpc to the source, 1 mas on the 
images corresponds to a projected separation of 10 AU. 
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6. VLBA IMAGE ANALYSIS 

The linear alignment of features in the VLBA images are 
interpreted as jets in an almost north-south orientation, 
in agreement with previous work (e.g. Mioduszewski et al. 
2001; Marti et al. 2001). 

The total flux densities recovered from the images track 
fairly well the flux densities measured with the VLA (Fig- 
ure 3). Perhaps a few hundred mJy is missing at 5 GHz (i.e. 
~ 5 %), although up to a few Jy (corresponding to several 
tens of percent) is not recovered at 22 GHz. However, it 
was difficult to measure accurately the flux densities of 
the different components from the images, mainly because 
of insufficient resolution. Hence it was difficult, especially 
for the images in which the r.m.s. noise is a significant 
fraction of the component flux density, to derive accurate 
spectral indices for the emission. At 22 GHz, the north- 
ern and southern components (initially unresolved) fade at 
similar rates, whereas at 5 GHz, the northern component 
fades more rapidly than the southern component. 

When considering the images, care should be taken not 
to overinterpret the details of the structures they show. 
In the 5-GHz and 15-GHz images, only the gross features 
should be relied upon, owing to the size of the scattering 
disk (16 mas at 5 GHz) or the high r.m.s. noise of the im- 
ages (126 mJy per beam in the 15-GHz image of September 
19). The fine structure in the 22-GHz images is certainly 
believable, owing to the low r.m.s. noise levels achieved 
and the high resolution of the images (a small beam size 
and a small scattering disk). Certainly on September 20, 
the shape of the jet, the extension to the north of the 
northern maximum, and the double peak to the south are 
believed to be real. 

6.1. Overall morphology 

Accepting that the extended emission we observed with 
the VLBA constitutes a jet, we now consider the nature 
of the jet; whether it consists of a single pair of ejecta, a 
continuous set of discrete ejecta, or a continuous hydrody- 
namic flow. 

The lower-frequency and hence lower-resolution (5-GHz 
and 15-GHz) images show the source to be dominated by 
two bright components. However, the highest-resolution 
image at 22 GHz reveals a sequence of several discrete 
knots, with the flux density being dominated by a sin- 
gle northern component (which we identify in §6.4.2 and 
§6.4.3 as the core with a northern extension) and a south- 
ern component which itself shows hints of being resolved 
into two distinct peaks. It seems most likely, therefore, 
that the jet consists of a series of discrete knots, but is 
dominated by two brighter components; one in the north 
and one in the south. In most of our subsequent anal- 
ysis and modelling we concentrate on only the brightest 
observed components seen at 5 and 15 GHz, since they ap- 
pear to dominate the flux density of the source and since 
we believe them to be uncontaminated by core emission 
(see §6.4.4). Figure 4 shows the movement of the bright- 
est emitting regions with time. 

Assuming the likely nature of the jet as a whole series 
of discrete ejecta, following the motions of the brightest 
components can help to determine some of the fundamen- 
tal jet parameters, e.g. the jet speed (§6.2). Measuring 
the evolution of the sizes of the brightest 22-GHz knots 



constrains the expansion speed of the plasmons and is evi- 
dence that we see discrete blobs. However, by considering 
the underlying nature of the jet, the positions of all the 
discrete knots seen in the 22-GHz image from September 
20 can be fitted to give information on the way the jet 
precesses (§6.4.3). 
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Fig. 3. — 22-GHz (top) and 5-GHz (bottom) VLBA component 
fluxes. The total flux detected with the VLBA tracks the flux 
detected with the VLA fairly accurately. Note how the northern 
component flux density drops dramatically at both frequencies be- 
tween September 19 and September 20, a decrease mirrored by the 
southern component at 22 GHz, but not at 5 GHz. The northern 
and southern components could not be distinguished at 22 GHz on 
September 18, and the flux density of the single measured compo- 
nent has been plotted for both, as it was not known how the flux 
density was split between them. 
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6.2. Jet speed 

The 15-GHz and 5-GHz data were phase-referenced, and 
analysis of the positions (depicted in Figure 4) of the two 
most prominent peaks in the VLBA images suggested that 
both the northern and southern components had non-zero 
proper motions. (We chose to ignore the 22-GHz data in 
these fits because these were not phase- referenced i.e. did 
not have absolute positional information and, for the rea- 
sons described in §6.4.4, because at this high frequency we 
believe there to be a significant core contribution.) We fit- 
ted a straight line to the measured positions at 5 GHz and 
15 GHz (both data sets taken together) as plotted in Fig- 
ure 4 (top and centre). We first constrained the gradient 
of the northern line to be zero (i.e. a horizontal line, indi- 
cating no systematic motion), which fit had a reduced x 2 
of 3.8. Fitting for an unconstrained gradient gave a signifi- 
cantly lower reduced x 2 of 0.64, with a best fit value for the 
proper motion of 2. 8 ±0.5 mas per day. This is smaller than 
the best-fit value for the southern component (8.4±1.7mas 
per day). Assuming the jet is two-sided (which we sub- 
stantiate in § 6.4), an analysis of the proper motions of the 
approaching and receding jets (e.g. Mirabel & Rodriguez 
1999) determines (3cosi = 0.50 ± 0.10 (where (3 = v/c 
and v is the intrinsic speed at which the knots move away 
from the core, and i is the inclination angle of the jet to 
the line of sight). This analysis also gives an upper-limit 
to the distance to the source of d = 35.7 ± 4.8 kpc and 
an ensuing upper limit to the inclination angle of the jet 
axis to the line of sight of i max = 59.8 ± 6.7°. Assuming 
a source distance of exactly 10 kpc, we can solve for the 
jet speed and inclination angle to get (3 = 0.56 ± 0.09 and 
i = 25.7 ±6.5°. 

A further check on derived proper motions is often made 
using a Doppler boosting analysis. Given the fading of the 
intensities of the plasmons with time (e.g. Figure 3 and 
§ 9.5) a simple Doppler boosting analysis is not valid. A 
refined application of this method (Mirabel & Rodriguez 
1999) which involves measuring the flux density in the 
approaching and receding knots at equal angular separa- 
tions from the core is in this case inapplicable because of 
the pronounced side-to-side asymmetry which is predicted 
from the jet speed and inclination angle derived above and 
which is seen in Figure 5. 

6.3. Plasmon size 

From the observations, it was possible to derive a maxi- 
mum size to the source components, which can be used 
to help estimate the magnetic field in the emitting re- 
gion. Since the images of Cygnus X-3 are highly scatter- 
broadened (Wilkinson et al. 1994), particularly at low fre- 
quencies, the apparent source size is dependent on fre- 
quency. Mioduszewski et al. (2001) reported that the ap- 
parent source size scales as 

= 448(^/1 GHz)" 2 09 mas. (1) 
which, within the errors, they found to be consistent with 
a v ^ scaling. Thus to determine the true source size, 
and to obtain the best resolution, we used our highest- 
frequency (22 GHz) images. The MPS task JMFIT was 
used to fit Gaussians to the two brightest knots in the im- 
ages, and the deconvolved sizes and orientations are shown 
in Table 4. The quoted sizes do not include a deconvolu- 
tion from the scattering disk, the size of which is given 
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Fig. 4. — Movement in Declination of the VLBA knots with time. 
Only Declination is plotted here because the jet axis is very close 
to north-south (§ 6.4.3 and Mioduszewski et al. 2001) and the knots 
are expected to move ballistically out from the core. Positions were 
measured by fitting one or more Gaussians to the emitting regions 
using the AIPS task JMFIT and finding the positions of the fitted 
peaks. Fitted positions are shown with crosses in Figure 2. For the 
points without a corresponding image in Figure 2, the position is 
taken from the obvious brightest point in a phase-referenced image 
which had not been self-calibrated. These images had not been in- 
cluded in Figure 2 as we were unable to produce good-quality images 
for these data sets. Uncertainties in the fitted positions are taken 
as the canonical value of 25% of the beamwidth. Uncertainties in 
the time are taken to span the course of each observation. Start of 
flare is taken as JD 2452167 ±0.2 (Trushkin et al. 2002), i.e. 2001 
September 14, but see §6.5 for discussion of when the ejection event 
occurred. 
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for the different observing frequencies in Table 3. The ob- 
served plasmon sizes increased with time, from 7.9^3Q to 
16. 9^3 6 mas; the errors on these sizes include the appar- 
ent elongation due to proper motion of the knots over the 
course of each observation. For the 8.5-hour observation of 
September 18, and a derived proper motion of 8.4masd _1 , 
the possible elongation of the knots is of order 3.0 mas. 
Thus the intrinsic size could be as small as 4.9 mas. What 
may complicate this simple measurement of the the plas- 
mon size however, is the possibility of there being multiple 
plasmons, as observed in the September 20 image, which 
sequentially make a contribution to the observed flux den- 
sity, but for which the VLBA has insufficient angular res- 
olution. Assuming that the measurements (incorporating 
the uncertainty due to the smearing effect of component 
proper motion) do reflect the plasmon expansion, the ex- 
pansion speed is then 4.6*° .1 mas per day, corresponding 
to a speed of O^Iq'^c for an assumed source distance 
of lOkpc. This is, reassuringly, fairly consistent with the 
plasmon expansion speed of 0.1 — 0.3c fitted by Atoyan & 
Aharonian (1999) for the 1994 outburst of GRS 1915+105. 

6.4. Is the jet one-sided or two sided? 

There is some discussion in the literature as to whether 
the observed jet in Cygnus X-3 is one-sided (owing to 
Doppler boosting effects reducing the flux density of the 
receding jet below the map noise as in Mioduszewski et al. 
(2001)) or two-sided (e.g. Marti et al. 2001). In the fol- 
lowing sections, we will try to assess whether the VLBA 
images of Figure 2 show evidence for a one-sided or a two- 
sided jet. 

6.4.1. Evidence from proper motions 

From the proper motion analysis of the 5-GHz and 15- 
GHz data in §6.2, the fits seem to favour the northern 
component seen at these frequencies being a moving jet 
rather than a stationary core. This argues for the ob- 
served VLBA jet being two-sided. This conclusion should 
be taken with the caveat that although the images were 
phase-referenced, they were also self-calibrated, a process 
which may in principle move component positions around 
by a fraction of a beam. Given that the beam size at 
5 GHz is between 13 and 16 mas and the northern compo- 
nent is seen to move by 12 mas over the course of the ob- 
servations, the observed movement could in principle have 
been caused by the self-calibration. We believe however 
that the systematic increase in declination of the northern 
component seen over a number of days at both frequencies 
is indicative of real proper motion. 

6.4.2. Core position 

Mioduszewski et al. (2001) showed convincing evidence 
that the northern component they observed during the 
1997 February flare was a core, rather than a receding 
counterjet to the extension seen to the south. The position 
of this core in their images appears to be 20 h 32 m 25 s .7733 
+40°57'27".965 (J 2000). Comparing this with the posi- 
tions in our images shows it to be almost coincident with 
the maximum in the northern component of the 22-GHz 
image from 2001 September 20. This would seem to argue 
that this component in our image (labelled "CI" in Fig- 
ures) is the core. However, Mioduszewski et al. stated 



that their phase referencing was not entirely successful, so 
there is some small uncertainty in their derived position 
(although their images prior to self-calibration showed that 
the northern component was stationary to within 3 mas). 
Moreover, the 1997 February outburst was 4.5 years ear- 
lier than that of 2001 September. If the system had a 
systematic velocity (owing to, for instance, a natal kick) 
of 100 km s -1 , it would move a distance « 95 AU in 4.5 
years. If this motion were in the plane of the sky, that 
would imply a shift in position of 10 mas at the distance 
of Cygnus X-3. That this should be along the jet axis is 
unlikely, but it should be borne in mind that the source 
core position might in principle have changed in the inter- 
vening time between the outbursts. 

6.4.3. Jet precession and modelling 

The final 22-GHz image (from September 20) shows a 
series of emitting knots, and the brightest northern com- 
ponent is coincident with the position of the core as deter- 
mined by Mioduszewski et al. (2001). We attempted to fit 
the precessing jet model of Hjellming & Johnston (1981) to 
this image, investigating the three most plausible core po- 
sitions, corresponding to the brightest point in the north- 
ern component (at 20 h 32 m 25 s . 77335 +40°57'27".9650), the 
flux density peak immediately south of that (at 20 h 32 m 25 s . 77315 
+40°57'27".9555), and the flux density peak at the centre 
of the jet flow (at 20 h 32 m 25 s . 77295 +40°57'27".9487), all 
given in J 2000 co-ordinates. These positions correspond 
to the points labelled "CI", "C2", and "C3" respectively 
in Figure 5. We used a minimisation routine based on the 
downhill simplex method (Press et al. 1992) in order to 
find the best agreement between the model and the im- 
age. The precessing jet model has eight independent pa- 
rameters: the inclination of the jet axis to the observer, 
i; the opening angle of the precession cone, ip; the angle 
through which the jet is rotated in the plane of the sky 
from an E-W alignment, \ (related to the position angle 
by P.A. = x + 90°); the source distance, d; the jet speed 
(3c; the precession period, P; the sense of precession s rot 
(clockwise or anticlockwise) ; and some initial phase of pre- 
cession, cj). We added in an extra parameter i v iew, the time 
after ejection at which the image was made. We measured 
the positions and flux densities at fifteen different points 
along the spine of the jet, and then weighted the points 
according to the corresponding flux densities. It was found 
that a simple minimisation of the sum of the weighted dis- 
tances between each measured jet position and the closest 
model point did not give reasonable solutions, since the 
minimised model parameters would predict radio emission 
from points well outside the observed jet. To improve the 
fits, additional terms were used in the function to be min- 
imised, penalising a fit if it predicted too many precession 
periods to have passed in the time between ejection and 
observation, and penalising it if it predicted radio emission 
outside a small box containing the observed jet. The min- 
imisation routine was run for almost a hundred thousand 
different sets of starting parameters, which were chosen 
in order to fully sample (at equal intervals) the available 
parameter space without taking a prohibitively long time 
to run. Results predicting (3 > 1, distances outside the 
range 5 < d < 15kpc (Predehl et al. 2000; Dickey 1983), 
or negative periods were disallowed. The sum of the dis- 
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tances from each jet position to the closest model point 
was computed for each of the converged results, and the 
results with the lowest such cumulative distances which 
also predicted a viewing time £ v iew < 8.5 days (the time 
between the detection of the initial rise of the radio flux 
density and the epoch at which the image was made) were 
then examined individually, to check that the results with 
the lowest cumulative distances did in fact produce good 
by-eye fits to the data. Upon examination, no satisfactory 
fits were found for core positions corresponding either to 
the centre of the jet flow (core position C3), or to the 
southern part of the northern component (C2). However, 
good by-eye fits to the data were found for a core position 
corresponding to the brightest point in the north compo- 
nent (CI). We therefore identify the brightest point in the 
northern component as the core of the system. 

For that northernmost core position, more good fits were 
found for clockwise rotation (s rot = —1) than for anti- 
clockwise rotation (s rot = +1), and the mean, standard 
deviation, and the maximum and minimum values for the 
satisfactory fits are given in Table 5. The sense of the jet 
precession can not be uniquely determined in this case, but 
the parameters ip, (3, P, and i v iew are similar for the two 
possible senses. <j) an d d are not uniquely determined, with 
<j> having a large standard deviation and the distance vary- 
ing between the limits allowed by the model-fitting code. 
Although the parameters i and \ appear very different in 
the two cases, we believe that this is due to the way the 
angles were defined in the model of Hjcllming & John- 
ston (1981). In changing the sense of jet precession from 
clockwise to anticlockwise, all angles will be measured in 
the opposite sense, which will affect the definitions of in- 
clination angle and position angle on the sky. Redefining 
*now = 180 — i and x ncw = X — 180 for the anticlockwise 
case gives i DCW = 10.4 ±3.8 and Xncw = 98.1 ±2.8, in good 
agreement with the parameters derived for the clockwise 
case. So, taking the core position to be at 20 h 32 m 25 s . 77335 
±40°57'27".9650, as allowed by the fits, we can constrain 
the system parameters to be ip = 2.4±1.2°, i = 10.5±4.2°, 
X = 98.9±2.9°, (3 = 0.63±0.09, and P = 5.34±0.67days. 
This gives a value /3cosz = 0.62 ± 0.11. Comparing this 
value with that of ftcosi = 0.50 ± 0.10 derived from the 
proper motion analysis of §6.2, we find a slight discrep- 
ancy, although the error bars do overlap. Our best model 
fits for both senses of precession are shown in Figure 5. We 
cannot robustly determine the sense of precession with this 
data. 

Thus we do see a two-sided jet on VLBA scales. How- 
ever, the jet speed [3 and the inclination angle i are such 
that the observed length of the northern jet is significantly 
shorter than that of the southern jet. The ratio of the arm 
lengths is 



1 + (3 cos i 
1 — (3 cos i ' 



(2) 



where / app is the length of the approaching jet and Z roc is 
the length of the receding jet, (3 is the jet-knot advance 
speed, and i the inclination of the jet axis to the line of 
sight. The fitted parameters in the precession model imply 
a length ratio of 4.3 ± 0.3. Figure 5 shows that the north 
jet is indeed shortened by a factor appropriate to these 
values of (3 and i, convincingly verifying the two-sidedness 
of the jet in Cygnus X-3 on this occasion. 



6.4.4. VLBA component spectra 

The flux densities of the source components (shown in 
the images of Figure 2) during the flare of 2001 Septem- 
ber were measured using the AIPS task JMFIT. The 5- 
GHz VLBA component lightcurves are shown in Figure 3. 
The total emission (the sum of the northern and southern 
components) was found to track the VLA flux density. At 
5 GHz, the integrated flux density of the northern compo- 
nent exceeded that of the southern component on Septem- 
ber 18 and 19, after which its flux density diminished 
more rapidly, and the southern component dominated. At 
22 GHz, the northern and southern components could not 
be distinguished on September 18. On September 19, the 
northern component dominated, and the flux densities of 
the two components were comparable on September 20. 
When the spectral indices of the VLBA components were 
determined, it was found that neither component had an 
flat or inverted spectrum between 5 GHz and 22 GHz. Fig- 
ure 5 shows that the core and the northern (receding) jet 
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Fig. 5. — Results of fitting the model of Hjcllming & Johnston 
(1981) to the 22-GHz VLBA image of September 20. The black dots 
overlaid on the VLBA image correspond to the derived positions of 
knots ejected every 0.1 days from the core. The bottom image is for 
clockwise rotation (s r ot = — 1) and the top image for anticlockwise 
rotation (s ro t = +1)- Both models fit all the major peaks in the 
southern jet, the direction of the northern jet, and the approximate 
lengths of the two jets. They make no allowance for the fading of 
jet knots with time. 
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are very close together owing to relativistic aberration, 
and therefore would not be resolved at 5 GHz. It is the 
jet which dominates the flux density at low frequencies 
(5 GHz and 15 GHz), but the core dominates at 22 GHz 
because it has a significantly less steep spectrum than the 
jet material. Therefore spectral indices between these fre- 
quencies will not (without significantly higher resolution) 
reveal the identity of the northern emission, be it core or 
jet. 

6.5. Is the ejection of the VLB A emitting regions 
coincident with the start of the VLA flare ? 

The time of the beginning of the radio outburst can 
be constrained to being on or prior to JD 2452167 (2001 
September 14), from the RATAN-600 data (Trushkin et 
al. 2002), and from Ryle Telescope monitoring (Pooley 
2001) 1 . The Ryle Telescope monitoring shows that the 
15-GHz flux density began to rise out of the pre-outburst 
"quench" period from 0.02 Jy on JD 2452162.5 to 0.1 Jy by 
JD 2452165. 

Extrapolating the positions of the VLBA knots at 15 GHz 
and 5 GHz (the absolute 22-GHz positions being inaccu- 
rate due to the lack of phase-referencing) in Figure 4 back- 
wards in time, we believe that those knots cannot have 
been ejected at the time when the low-resolution tele- 
scopes (the RATAN-600 and the VLA) detected the radio 
flux density beginning to rise. Extrapolating to a precise 
time of ejection is difficult, owing to the large uncertain- 
ties on the measurements, but it seems that the northern 
and southern components would have been coincident per- 
haps 2.5±0.5 days after the RATAN-600 saw the radio flux 
density start to rise at the beginning of the outburst. This 
is 1.3 ± 0.5 days prior to the first VLBA observation on 
2001 September 18. By way of explanation, we propose 
a scenario in which the core flared and only expelled the 
jets later. This we explore in a forthcoming paper. Note 
that this is the opposite order of events to that proposed 
by Atoyan & Aharonian (1997), who suggested a scenario 
whereby the radio flux density flared well after the ejec- 
tion of radio clouds, once they had expanded sufficiently to 
become transparent at radio frequencies. Alternatively, it 
may be that the jet speed increased over the course of the 
first few days, although we note that once VLBA monitor- 
ing began (on September 18) the positions of the knots are 
consistent with a constant separation velocity, suggesting 
that their motion has always been ballistic. 

6.6. Comparison with previous work 

The most detailed previous VLBI images of Cygnus X- 
3 were made by Mioduszewski et al. (2001). They made 
three observations at 15.365 GHz with the VLBA at epochs 
2, 4, and 7 days after the peak of a 10 Jy radio flare in 1997 
February. Mioduszewski et al. detected what they inter- 
preted as a stationary (to within 3 mas) core component 
and a continuous jet extending to the south, reaching a 
maximum observable extension of 120 mas four days after 
the flare, which faded beyond the detection limits by the 
last observation. They took the time of ejection of the 
radio-emitting plasma to be the start of the outburst as 
detected by the Ryle Telescope and the Green Bank In- 
terferometer (GBI), and from the measured lengths of the 

1 http://www.mrao. cam. ac.uk/~guy 



jet at different epochs, derived proper motions, which con- 
strained the jet speed to be /3 m i n > (0.78 — 0.75), and a 
weak upper limit to the inclination to the line of sight to be 
i < 78 . Then, from the ratio of the integrated flux den- 
sity in the southern jet to the integrated noise in the north 
taken over a smaller area (to account for slower apparent 
motion of the receding jet), they derived a constraint on 
,Scosz > 0.806. We believe that there is some scope for 
error in their proper motion analysis, since it rested on 
the length of the jet, rather than on tracking individual 
well-identified knots of radio emission. The length of a 
jet that may be imaged will be affected by the map qual- 
ity, which was poor during at least the first epoch of their 
observations, owing to bad weather in the southwestern 
USA. Moreover, it is not clear that the radio plasma was 
ejected at the start of the Rylc/GBI flare (see §6.5). Un- 
certainty in the ejection time would then lead to an error 
in the proper motion and the derived jet speeds. 

A comparison may also be made between the precession 
parameters derived for the 1997 flare and those we derived 
in §6.4.3. Their best fits yielded a precession period > 60 
days, a jet speed (5 > 0.81, a precession opening angle of 
<12°, an inclination angle of <14° to the line of sight, and 
anticlockwise precession. The inclination angle, jet speed 
and opening angle derived by Mioduszewski et al. (2001) 
were similar to those we derived. 

Because the derived precession period of the 2001 event 
is comparable to the monitoring timescale (c.f. the 1997 
event, when the monitoring sampled less than one tenth 
of the best-fit precession period), we have confidence in the 
limits we have placed on this parameter. However, the fits 
of Mioduszewski et al. robustly require a long precession 
period (> 30 days), whereas ours is of order ~ 5.3 days. 
It is thus hard to avoid the conclusion that the precession 
period has changed by a significant factor between these 
two events. 

From the above, we believe that the jet speed we derive 
from proper motion measurements is consistent with that 
estimated by Mioduszewski et al. (2001). In both sets of 
observations, the southern jet appears to be brighter and 
faster-moving than any northern jet and in both sets of 
data there is evidence for precession (albeit with rather 
different precession parameters). In addition, the nature 
of the jet itself appears to be different, in the sense that the 
2001 September jet seemed to be composed of a few identi- 
fiable knots, whereas the 1997 February jet appeared more 
like a continuous stream of radio plasma. On the other 
hand, intrinsic differences could, in principle, account for 
the difference in the observed nature of the jet. The out- 
burst of 1997 February was shorter in duration than that 
of 2001 September. The integrated flux density died away 
more quickly in the 1997 flare (falling to below 1 Jy at 
15 GHz ~ 3 days after the peak of the flare, as compared 
to ~ 5 days afterwards in 2001), after which it was still 
possible to make VLBA observations. The dynamic range- 
limited images of Mioduszewski et al. therefore have much 
lower r.m.s. noise levels than those in Figure 2, so lower- 
level continuous emission (as opposed to the discrete knots 
we observed) could be detected. This could account for our 
non-detection of continuous emission, although the com- 
parison of VLA and VLBA flux densities suggest that we 
are not missing significant components in our VLBA im- 
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ages, at least at low frequencies. 

A comparison may also be made to the flare of 2000 
September 15 observed by Marti et al. (2001). They ob- 
served Cygnus X-3 with the VLA in its most extended 
"A" configuration 36, 51, and 66 days after the start of a 
radio outburst detected with the GBI. By this latest time, 
the source had faded back to quiescence. They imaged a 
two-sided ejection on arcsecond scales, with a dominant 
core, a prominent northern ejection and a strong hint of 
extension to the south. The flare observed with the GBI 
was less powerful than either the 1997 or the 2001 flare, 
peaking at only ~ 7 Jy at 8.3 GHz. 

Marti et al. (2001) derived proper motions by assum- 
ing the ejection date to be the start of the flare observed 
with the GBI, and obtained a value of [3cosi = 0.14 ± 
0.03. In fact, their plot of the GBI data shows several 
flares of amplitude > 4 Jy over the 20 days following their 
assumed ejection date (after which no data is plotted, 
since GBI monitoring of Cygnus X-3 appeared to cease 
on JD 2451821). This opens up two possibilities. It is not 
clear which (if any) of the plotted flares was associated 
with the actual ejection event corresponding to the jets 
they observed. If the last one (JD 2451822) was the flare 
responsible, their third image showing the arcsecond-scale 
jets would have been taken only 46 days after the flare. 
This would lead to increased values for the proper mo- 
tion of the two jet components, and a derived jet speed of 
(3 = 0.68. It is also possible that the ejection occurred af- 
ter their GBI data stops, in which case the derived proper 
motions and speeds would be even greater. A speed of 
(3 = 0.68 is much more consistent with that derived in 
this paper. It is not clear however that the expansion rate 
would be uniform and it is possible that the jet would slow 
down as it impinged on its environment. While this effect 
might be negligible over a few days, it might be signifi- 
cant on a timescale of months. However, their value of 
(3 cos i = 0.12 ± 0.01 based on the flux density ratio of the 
2 jets interpolated to equal angular distances from the core 
would seem to be more robust than the value derived from 
their proper motion analysis. 

Marti et al. explained the one-sided VLBA jet observed 
by Mioduszewski et al. (2001) by invoking jet obscuration 
by an absorbing medium on VLBA scales, in the form of 
a free-free absorbing disk-like wind. They suggested that 
the northern jet was the approaching jet, and the south- 
ern jet was receding from us. For the 15-GHz images of 
Mioduszewski et al., they derived an optical depth at 50- 
lOOmas from the core of — 1.15 ± 0.14. Since their 
quoted free- free absorption coefficient scales as i/~ 2,08 , that 
would imply an optical depth of 0.53±0.08 at 22 GHz. This 
would be insufficient to prevent us from seeing a bright ap- 
proaching northern jet in our images. Since our observa- 
tions strongly suggest a receding northern jet, it would be 
hard to reconcile their obscuration explanation with our 
data. 

We made further observations of Cygnus X-3 134 days 
after the detected start of the 2001 September flare, on 
2002 January 25, whose detailed analysis we will present 
in a forthcoming paper. As with the observations of Marti 
et al., the VLA was in "A" configuration, and the source 
had faded back to a quiescent flux density of « 100 mJy 
at 43 GHz, although it displayed significant variability on 



a timescale of ~ 10 minutes. We subtracted out the time 
variability, and then imaged the core. Fitting that with a 
Gaussian and subtracting this from the uw-data set, there 
was no evidence of any extension similar to that seen by 
Marti et al., to an r.m.s. limit of 0.8 mJy per beam. We 
note, however, that this is as expected, since Marti et al. 
(2001) detected northern and southern extensions 66 days 
after the assumed ejection event (see discussion above) at 
levels of 1.45 mJy and 0.89 mJy respectively. Assuming 
that the flux densities of these extensions decreased with 
time, then 134 days after the start of the outburst, such 
extensions would be undetectable at our r.m.s. limit. 

7. RADIO SPECTRA 

7.1. General characteristics 

The high-frequency parts of the spectra shown in Fig- 
ure 6 were fitted with single power laws and the resultant 
spectral indices are shown in Figure 7. The spectral in- 
dex evolved during the course of the flare, from a value of 
0.44 ± 0.06 on September 18 to a value of 0.61 ± 0.08 by 
September 23. Traditional models of synchrotron radia- 
tion cannot explain such a spectral steepening, although 
this phenomenon has been seen before, in the giant radio 
outburst of 1972 (Aller & Dent 1972; Hjellming & Balick 
1972; Peterson 1973), and in subsequent outbursts (e.g. 
Seaquist et al. 1974; Geldzahler et al. 1983). 

The spectra also show that by September 24 (if not 
earlier), the spectra began to deviate from pure power 
laws, and rapid variability began to affect the source flux. 
Two different scans were made at 8.4 GHz on Septem- 
ber 24, each of duration 2 minutes, and separated by 30 
minutes (this was the only epoch when separated scans 
were made) ; in this time, the source flux density increased 
from 618.2 ± 1.2 mJy to 656.5 ± 1.2 mJy. The main out- 
burst dominated the source flux during the first few days 
of observation, but as the outburst decayed, the source 
behaviour was affected more by core variability or minor 
ejections. 

7.2. What causes the low- frequency turnover? 
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Fig. 6. — Evolution of the radio spectrum of Cygnus X-3 
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The low-frequency (1.4-GHz to 74-MHz) measurements 
indicate a turnover in the radio spectrum. We now con- 
sider various possible mechanisms for creating such a turnover, 
and assess the plausibility of each one. Such low-frequency 
turnovers in other sources are generally ascribed to one of 
four distinct mechanisms (Kellermann 1966; Gregory & 
Seaquist 1974). These are: synchrotron self-absorption 
(SSA), a low-energy cutoff in the electron energy distri- 
bution spectrum, the Tsytovich-Razin (T-R) effect, and 
free-free absorption (FFA) by thermal plasma, either with 
the plasma uniformly mixed with the radiating electrons 
or as intervening material along the line of sight to the 
source. We now discuss these four possibilities in turn. 

7.2.1. Synchrotron self- absorption 

Well below the turnover frequency z/ m , where the op- 
tical depth t 3> 1, a synchrotron self-absorbed spectrum 
scales as S v oc v b / 2 (e.g. Rybicki & Lightman 1979). In 
the optically-thin region however, where r <C 1, the scal- 
ing is S v oc v~ a (a > 0). The width of the turnover 
region depends on the structure of the source. Unfortu- 
nately, we do not sample the low-frequency spectrum well 
enough to determine whether the terminal low-frequency 
spectrum scales as v b / 2 . However, we clearly observed that 
the turnover moved to lower frequencies with time. If the 
synchrotron self- absorption paradigm is correct, this would 
be consistent with the observed expansion of the plasmons, 
which would cause the magnetic field within the plasmons 
to decline, and thus lower the turnover frequency. Mag- 
netic flux freezing implies that B/p oc Al (Landau & Lif- 
shitz 1963), where p is the plasma mass density, B is the 
magnetic flux density, and Al is the length of a fluid el- 
ement along the magnetic field line. Assuming that the 
mass of the plasma cloud is constant and that there is no 
strong turbulence, Al(t) oc R(t), and B oc R~ 2 . However, 
if turbulence develops in the plasma, the fluid elements are 
stretched by turbulent eddies in addition to the expansion 
of the cloud, so Al oc R 1+a , where a in this case is a fac- 
tor describing the degree of turbulence, typically taken as 
of order 1 for a turbulent plasmon with a correspondingly 
tangled magnetic field. In this case, B oc R~ 2+a « R~ x . 
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Fig. 7. — Evolution of the spectral index of Cygnus X-3. The 
spectral index was fitted with a simple power-law for the portion 
of the spectrum that was well above the low-frequency turnover 
(> 8.4 GHz for September 18 and 19, > 5 GHz for September 21 
and 22, and > 1.4 GHz for September 23 and 24). 



So as the plasmon expands, the magnetic field declines, 
and the turnover frequency decreases. Thus the SSA sce- 
nario is consistent with our data. 

7.2.2. Low-energy cutoff in the electron spectrum 

We modelled the electron distribution as a power law. 
If that power law were to extend only down to some mini- 
mum energy, below which the electron distribution turned 
over, producing a cutoff in the electron distribution func- 
tion, that would produce a low-frequency turnover in the 
radio spectrum. This turnover would move to lower fre- 
quencies as the energies of the individual electrons de- 
creased due to adiabatic expansion of the plasmon, moving 
the cutoff in electron energies to lower energy. We cannot 
rule out this possibility with our data. Depending on the 
abruptness of the cutoff, observing a v 1 / 3, spectrum below 
the turnover would identify this scenario. 

7.2.3. The Tsytovich-Razin effect 

The refractive index, n, of a plasma containing free elec- 
trons and ions may differ from unity and signals then travel 
at the phase velocity c/n. In a medium where the refrac- 
tive index is less than one, beaming from relativistic elec- 
trons is suppressed, as the angle into which the radiation 
is beamed is given by 

l/ 7 = ^l-n 2 [3 2 . (3) 

At low frequencies [y <C v p , where v v is the plasma fre- 
quency), this effect is dominant, and the synchrotron spec- 
trum is cut off due to beaming suppression. At higher 
frequencies, this effect is unimportant, since 

rf-i-e)*. <4, 

The frequency in GHz at which this effect becomes signif- 
icant is given by 

i/~2x 1Q- I4 n th /B, (5) 

where B is the magnetic field in Gauss and nth is the 
density of ambient thermal matter in m~ 3 . In §8.1 we 
estimate an upper limit on the magnetic field to be B ~ 
0.15 G; thus if the T-R effect were the mechanism by which 
the spectrum turns over, a rather high (upper limit to 
the) density of 8 x 10 12 m~ 3 of thermal plasma would be 
required for a turnover frequency of 1 GHz. At these densi- 
ties, free-free absorption would also be expected to become 
important. 

7.2.4. Free-free absorption 

We consider two possible geometries for free- free absorp- 
tion. The absorbing thermal material could either form a 
screen enclosing but not permeating the source, or it could 
be mixed in with the synchrotron-emitting plasma. 

In the former case, if the plasmons were expanding into 
a region of thermal material, then the flux of the receding 
component would fall off more rapidly with time than the 
flux of the approaching component, since the length of the 
absorbing screen would increase as the receding plasmon 
moved away from the core. The length of the absorbing 
screen between the observer and the approaching compo- 
nent would decrease with time, so the turnover would move 
to lower frequencies with time as this component began to 
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dominate the observed flux density. In fact, we see ev- 
idence that the 5-GHz flux density of the northern (jet) 
component falls off rapidly with time after September 19, 
whereas that of the southern component does not diminish 
so rapidly. This is then in principle a plausible scenario. 

For thermal material to be mixed with the synchrotron- 
emitting plasma, the material would either need to be 
present at the time the plasmon was ejected, or would 
have to be entrained during the subsequent motion of the 
plasmon through the surrounding medium. In order to 
create a turnover at frequencies as high as ~ 3 GHz by 
September 18, the initial rate of entrainmcnt would have 
to be extremely high. To cause the observed decrease in 
turnover frequency, the rate would then have to be quite 
accurately tuned. We do not favour this explanation. 

7.2.5. Discussion on the origin of the turnover 

The change in turnover frequency with time places more 
stringent constraints on the mechanism responsible for the 
turnover. For the T-R effect, the change in turnover fre- 
quency with time would require a change in n t h/B. For 
FFA, it requires a decreasing density with time, either 
due to the source moving into lower-density regions as it 
moves outwards, or due to the expansion of thermal mate- 
rial mixed in with the synchrotron-emitting electrons. For 
SSA, expansion and a corresponding decrease in the mag- 
netic field lead to a decrease in the absorption coefficient 
and therefore in the turnover frequency. Multi-frequency 
imaging with sufficiently high resolution at or below the 
turnover frequency could in principle distinguish defini- 
tively between the mechanisms. 

7.3. High-frequency spectral evolution 

In order to explain the high-frequency spectral evolu- 
tion, in § 9, we develop a model which factors in light travel 
time effects both within and between components, which 
we believe explains the basic characteristics of the spectral 
evolution, e.g. the observed timescale. But first we derive 
some basic parameters of the source in § 8. 

8. CONSTRAINING THE MAGNETIC FIELD AND LORENTZ 

FACTORS 

8.1. Magnetic Field 

If self-absorption were the mechanism responsible for 
the turnover, precise size constraints from VLBA measure- 
ments would allow the direct determination of the mag- 
netic field. If, however, another process were responsible, 
the self-absorption turnover frequency v m would be at a 
lower frequency than the observed turnover in the spec- 
trum. Thus an extrapolation of the power-law spectrum 
back to the predicted self-absorption turnover frequency 
would give a higher extrapolated turnover flux density S m 
than is observed. Since B cx S 1 " 2 ^, assuming that the ob- 
served turnover frequency and flux density are due to SSA 
still gives an upper limit to the magnetic field strength. 

Assuming that the turnover in the spectrum of Cygnus 
X-3 is due to SSA, then given the frequency, f m , and flux 
density, S m , of the turnover, and the angular size 9 of the 
emitting region, synchrotron theory gives an expression for 
the magnetic field in that region (Marscher 1983). 

B = lO- 5 b(a)0 4 v 5 m S m 2 G, (6) 



where S m is measured in Jy, 9 in milliarcseconds and v m 
in GHz. The parameter b(a) is a slowly- varying function 
of the high-frequency spectral index with a value of w 3.4 
for a = 0.6, and S is the beaming parameter 



r(l -/3cos*)' v ' 

where T is the bulk Lorentz factor, (1 - /3 2 ) -1 / 2 . The 
rcdshift z is zero. 

Using the VLBA observations to constrain the number 
and sizes of synchrotron-emitting components, it is possi- 
ble to fit spectra to those components and derive upper 
limits to their magnetic fields. Unfortunately, the strong 
dependence of the magnetic field on the angular diameter 9 
and the turnover frequency v m means that the uncertainty 
in the derived field is quite high. It also makes it extremely 
important to define correctly the angular diameter 9. The 
radiative transfer equation used to derive Equation 6 as- 
sumes spherical components. The fitting task in AIPS 
used to derive the sizes models the components as ellipti- 
cal Gaussians. Marscher (1987) gives an expression for the 
value 9 S of a sphere which subtends the same solid angle 
as an ellipse of major and minor FWHM axes #c a an d $G b 
as 

9 S = 1.8y/e G jG b - (8) 

These equivalent angular sizes are listed in Table 4. Us- 
ing Equations 6 and 8, together with the fitted turnover 
frequency of 3.01 ± 0.06 GHz, the turnover flux density of 
15.1 ± 0.5 Jy, the spectral index of 0.4, and the measured 
core size of 7.88 ± 0.51 mas on September 18 (when the 
VLBA observations suggest that the source consisted of 
a single component; we note that even if it did not, the 
individual component sizes would then be smaller than 
the fitted size, which we could therefore still use to place 
a reasonably firm upper limit on the magnetic field), an 
upper limit on the magnetic field of 1.5 ± 0.6 x 10 _1 G 
(1.5 ± 0.6 x 10~ 5 T) is derived. 

This is our best estimate of the magnetic field strength. 
It could, however, be at least an order of magnitude lower. 
The proper motion of the plasmons during the course of a 
VLBA observation would elongate the observed knots in 
the images and cause us to overestimate their sizes. If this 
effect were taken into account, the derived magnetic field 
might be at least as low as 1.94 ± 0.96 x 10~ 2 G. 

8.1.1. Equipartition? 

If the particle and magnetic field energies in the source 
components were close to equipartition, it would be pos- 
sible to estimate the minimum-energy magnetic field in 
the plasmon. Following the derivation set out in Longair 
(1994), and using measurements of the flux densities and 
plasmon sizes (Table 4) from September 18, and assuming 
a factor of 100 times more energy in relativistic protons 
than in electrons, this assumption gives a minimum-energy 
magnetic field of 6 x 10~ 2 G, corresponding to a minimum 
energy in the plasmon of 3 x 10 34 J (an energy density of 
4 x 10 -5 Jm -3 ). Since our upper limit on the magnetic 
field is greater than the minimum-energy field, our mea- 
surements are not necessarily inconsistent with equiparti- 
tion (although we note that there is no compelling reason 
why the equipartition assumption should hold in such a 
manifestly non-steady state system). 
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8.1.2. Comparison to other magnetic field estimates 

Most derivations of the magnetic field in microquasar 
jets come from equipartition arguments, despite the uncer- 
tainty as to whether equipartition is valid in these systems. 
Spencer ct al. (1986) used equipartition arguments to put 
a lower limit of 0.08 G on the magnetic field of the emit- 
ting components during the 1983 September outburst of 
Cygnus X-3. Equipartition assumptions applied to other 
sources have also led to minimum-energy magnetic fields 
of order 0.01 - 1 G (for instance GRS 1915+105 (Atoyan 
& Aharonian 1999), LS5039 (Paredes et al. 2000), and 
Scorpius X-l (Fomalont, Gcldzahler & Bradshaw 2001)). 

A few magnetic field estimates are available in the lit- 
erature which do not rely on equipartition. Ogley et al. 
(1998) used electron ageing arguments to put an upper 
limit on the magnetic field in the jet of Cygnus X-3 in qui- 
escence of B < 20 G at 5 x 10 12 cm from the core. Paredes 
et al. (2002) considered inverse Compton processes from 
a conical jet to constrain the magnetic field at 1000 AU 
from the core of LS 5039 to lie between 0.06 and 128 G. 
Thus our estimate of the magnetic field in Cygnus X-3 is 
consistent with derived values in both this system and in 
other microquasars. 

8.2. Lorentz factor 

To a good approximation the radiation observed at a 
given frequency v probes electrons of Lorentz factor 

{2irm c v\ 1/2 

where B is the magnetic field in T, v the observing fre- 
quency in Hz, and the constants arc in S.I. units. In the 
case of adiabatic expansion, the magnetic field decreases as 
the source expands, and observations at a fixed frequency 
sample electrons with higher Lorentz factors. 

From the derived upper limit to the magnetic field on 
September 18 of 0.15 G, and the maximum and minimum 
observing frequencies on that day (43.3 and 1.4 GHz re- 
spectively), the lower limits to the Lorentz factors of the 
electrons probed by the measurements are 321 and 58 re- 
spectively. 

8.3. The B — 7 plane 

Synchrotron theory, taken together with our data, can 
be used to derive constraints on the area of the B — 7 plane 
(Figure 8) occupied by the source in 2001 September. 

8.3.1. Synchrotron self-absorption turnover 

The upper limit on the magnetic field derived in §8.1 
from synchrotron self-absorption arguments is shown as 
line 1 in Figure 8. 

8.3.2. Energetics 

If we assume that the energy ejected during the out- 
burst cannot greatly exceed the energy emitted by a body 
radiating at the Eddington luminosity over the rise time 
of the outburst (taken to be of order 1 day), then an up- 
per limit on the total energy in the plasmon should be 
of order ~ 10 36 (M/M Q ) J. This energy is contained in 



the relativistic particles and the magnetic field. An up- 
per limit to the magnetic field can be found if all the en- 
ergy is assumed to be in the field rather than the parti- 
cles. Knowing the volume of the plasmon (assumed spher- 
ical) on September 18, this leads to a limit on the en- 
ergy density of 10 _3 (M/Mq) Jm -3 . By measuring the 
Dopplcr shifts of X-ray spectral lines, Stark & Saia (2003) 
placed an upper limit on the mass of the compact object 
of 3.6 Mq (subject to certain assumptions). Thus an upper 
limit of 6xlO- 1 (M/M Q ) 1 ^ 2 = 1.1 G may be placed on the 
magnetic field (shown as line 2 in Figure 8). 

8.3.3. Absence of catastrophic synchrotron losses 

The fact that the spectrum shows no noticeable break 
up to 100 GHz implies that synchrotron losses cannot be 
important up until at least the time of the last observation 
on September 24. The flare began on September 14. Thus 
the synchrotron loss timescale i syn must exceed ten days. 
This puts a limit on B 2 -f, namely 

B 2 7 < 8.9 x 10 2 G 2 . (10) 

This is shown as line 3 in Figure 8. 

8.3.4. Lack of leakage 

Another constraint may be derived by considering the 
synchrotron-emitting electrons. To emit synchrotron ra- 
diation, they must be gyrating in a magnetic field. In or- 
der for the electrons to remain within the plasmon (where 
the magnetic field is high enough for them to produce de- 
tectable radio emission), their gyroradius must be signifi- 
cantly smaller than the size of the plasmon. If it becomes 
too large, electrons will begin to leak out of the plasmon 
and will no longer produce detectable emission. The gy- 
roradius r g of an electron of Lorentz factor 7 in a field B 

Requiring this to be less than some fraction / of the size 
of the plasmon gives 

7m c 



eB 



<f0d, 



(12) 



where 9 is the observed angular size of the plasmon and d 
is the distance to the source (taken as lOkpc). This limit 
is shown as line 4 in Figure 8. 

8.3.5. Region probed by the measurements 

Since the electrons radiate at particular frequencies, Equa- 
tion 9 can be used to find the exact lines on the B — 7 plane 
probed by the observations at each frequency (shown as 
lines 5 in Figure 8). 

8.4. Number and energy densities 

From Longair (1994) and integrating over a power-law 
distribution of electrons with exponent p ^ 2 (in acknowl- 
edgement of the terminal high frequency spectral indices 
we observe, even though p = 2 is the most frequently- 
quoted case as it simplifies the mathematics) an expres- 
sion for the number density of emitting electrons may be 
derived. Modelling the emitting region as a sphere of uni- 
form density, of radius r = j9d, the number density (in 
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m 3 ) of emitting electrons is 

2.344 x 10 - 25 (p - l)d0 3 a{p) (w) ^ ( 1 ' 253xl ° 37 ) ^ 

(13) 

where d is the distance to the source in metres, is the 
angular size of the source in radians, S v is the observed 
flux density (measured in Wm _2 Hz _1 ) of the source at 
frequency v (Hz), B is the magnetic field of the source in 
G, m c is the mass of the electron in kg, c is the speed 
of light (ms -1 ) and a(p) is a slowly-varying function of p 
which is tabulated in Longair (1994). 7 m i n and 7 max are 
the minimum and maximum Lorentz factors over which 
the electron spectrum is known to extend, and are thus 
the limits of integration used in calculating Equation 13. 

From the September 18 data, Equation 13 gives a lower 
limit on the electron density of 9.2 x 10 7 m~ 3 , although we 
note that this is uncertain for many reasons not least the 
limits of integration used (we used only the observed range 
of Lorentz factors listed in §8.2) and the uncertainty in the 
sizes of the plasmons and hence the derived magnetic field. 
A lower limit on the energy density in relativistic electrons 
therefore is 9.8 x 10~ 4 Jm~ 3 . An upper-limit on the energy 
density in the magnetic field is estimated (from the upper 
limit to the field derived in §8.1) to be 9.0 x 10~ 5 Jm~ 3 
on September 18. This gives a ratio of energy in particles 
to energy in the magnetic field of > 10 for that day. The 
magnetic field would then be a factor 3 less than predicted 



by equipartition arguments. 

9. HIGH-FREQUENCY SPECTRAL EVOLUTION 

The evolution of synchrotron-emitting plasma will be 
influenced by adiabatic expansion of the emitting region 
(van der Laan 1966; Scheuer & Williams 1968; Hjellming 
& Johnston 1988). The energy loss rate for a gas of rela- 
tivistic particles in an expanding volume is 

£=-±(V-v)£, (14) 

where E is the particle energy and v is the particle velocity. 
For a uniformly-expanding sphere of radius R, the radial 
velocity distribution in the sphere is 

v = R{r/R), (15) 

so that 



The energy loss rate is proportional to the energy of the 
particle, so that the expansion is self-similar, and the spec- 
tral shape is therefore preserved. 

We observed an initial high-frequency spectrum S{v) oc 
^-o.43±o.o6 w j 1 i c } 1 steepened between September 18 and 
September 21 to approximately S(v) oc ^~ ' 6 , and then 
remained consistent with this terminal value for the re- 
maining three days of observation. To cause the spectrum 
to steepen, some form of energy-dependent loss mechanism 
must be invoked. We now consider several possibilities for 
this. 



CO 




Fig. 8. — Constraints on the B — 7 plane. 1 is the upper limit 
on the magnetic field from SSA (Equation 6) and the size constraint 
of 7.7mas from VLBA observations (§8.3.1). 2 is the constraint 
from energetics considerations assuming a 3.6 solar mass compact 
object (§ 8.3.2). 3 is the synchrotron ageing limit from Equation 10 
(§ 8.3.3). 4 is the limit from the requirement for the gyroradius to 
be less than 10% of the plasmon size, for an angular size of 7 mas 
(§ 8.3.4). 5 are the loci in the plane probed by observations at three 
given frequencies (Equation 9). Dashed lines are for 99 GHz, dot- 
dashed lines for 8.4 GHz and dotted lines for 330 MHz. The arrows 
show which side of the constraint lines are valid parameter space. 
The shaded areas denote parameter space which is excluded. 



9.1. Synchrotron losses 
The energy loss rate for a synchrotron-emitting electron 



is 



E = --CT T ct/ ma g7 2 - 



(17) 



Therefore the synchrotron loss timescale is given by 

^ 55 1.3 x l0-»ffl 7 SeC ° ndS ' (18) 

where B is measured in Gauss. The upper limit on the 
magnetic field of order 10 _1 G (§8.1) leads to a lower limit 
on the synchrotron loss timescale t syn of order 11 years. So 
synchrotron losses cannot be responsible for the steepen- 
ing over a timescale of order days. Moreover, synchrotron 
losses are predicted to create a break in the spectrum 
which moves to lower frequencies with time. We do not ob- 
serve a break in the high-frequency spectrum, but rather a 
smooth spectrum with a single spectral index all the way 
up to 100 GHz, albeit one which initially varies with time. 

9.2. Bremsstrahlung losses 

The bremsstrahlung energy loss rate for an ultrarela- 
tivistic electron in a fully-ionised hydrogen plasma of num- 
ber density rin m -3 is 

E = -7 x 10" 23 n H (log7 + 0.36)S (19) 

which leads to an energy-loss timescale of 

t brom = 6.1 x lO 21 ?^ 1 seconds (20) 

for particles of Lorentz factor 10 2 . In order to give a cool- 
ing time of order days, a proton density of 1.6 x 10 16 m -3 is 
required. This is nine orders of magnitude greater than the 
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electron number density inferred in §8.4 and would imply 
a total baryonic mass almost equal to a solar mass in the 
plasmon. This proton density would give rise to a very 
large emission measure, and hence the free-free opacity 
would be very high, and we would not see the observed syn- 
chrotron emission. We therefore rule out brcmsstrahlung 
as the loss mechanism responsible for the spectral evolu- 
tion. 

9.3. Inverse Compton radiation losses 

An energetic electron in a radiation field of energy den- 
sity [/rad loses energy at the rate 

E = -§<7 T c£/ rad7 2 . (21) 

Hence the cooling time for inverse Compton radiation is 

7 _. on „,„6rr-l..-l (22) 



he = ~i^ 3.1 x 10 6 U-^-\ 



1 

For a particle of Lorentz factor 7 ~ 10 2 , a cooling time of 
order a few days would require a radiation field of order 
C^rad ~ 10 _1 Jm~ 3 in soft photons. Given this energy den- 
sity, a source of photons local to the plasmons would be 
required, since background radiation fields (the CMB, the 
interstellar background etc.) would not have the required 
intensity. 

If the electrons inverse-Compton scatter off an external 
source of photons of luminosity L emitted a distance r 
away, then if the plasmon is not moving relativistically, 

UnA = d^' (23) 

where all quantities given are in S.I. units. Inserting the 
energy density we require in order to have an inverse Comp- 
ton cooling time of order a few days gives a constraint on 
the photon luminosity we need and the distance of the 
plasmon from that source, 



^Edd 



6 - 6x1 ° 2 (kkJau) 2 ' ™ 



where LEdd is the Eddington luminosity for a ten solar 
mass companion star. Thus, we do not consider this as 
the likely dominant loss mechanism, since the companion 
star's luminosity would have to exceed its Eddington lu- 
minosity by at least two orders of magnitude. 

Fender et al. (1997) reported that inverse Compton losses 
were initially dominant over expansion losses at 15GHz 
during a flare in 1994 February. Their model showed that 
the relative significance of the inverse Compton losses then 
declined to ~ 10% of the expansion loss level by three days 
after injection. We note that they assumed a significantly 
slower jet velocity (~ 0.3c) than we found in § 6.4.3. They 
also proposed that radiation from the stellar wind of the 
companion star would contribute to the radiation field (a 
plausible assumption if the companion is indeed a Wolf- 
Rayet star). The combination of these two effects could 
in principle account for the lower significance that we at- 
tribute to inverse Compton losses. 

9.4. Leakage 

Energy-dependent escape of relativistic electrons from 
the plasmons was considered by Atoyan & Aharonian (1999) 
in order to explain the spectral steepening observed in the 



1994 outburst of GRS 1915+105. They also put forward a 
second possibility, that the steepening was caused by the 
continuous injection of electrons with a gradually steep- 
ening electron energy spectrum. Our VLBA observations 
of Cygnus X-3 seem to indicate a discrete set of ejecta, so 
we discount the continuous ejection scenario on observa- 
tional grounds. We now consider the first possibility, that 
of energy-dependent escape. 

Electrons within the plasmon will be scattered by tur- 
bulent magnetic fields and diffuse through the cloud as a 
result. This will cause diffusive leakage from the surface 
of the plasmon on a timescale R 2 /D where R is the plas- 
mon size and D the energetic electron diffusion coefficient 
which is D = Ac/3 where A is the mean free path. In 
quasilinear theory the smallest mean free path allowed is 
the particle gyroradius, which corresponds to a turbulent 
field of order the mean field, 



7mc 
~~e~B 



1.7 x 10~ 3 ^-, 
B 



and hence 



1.4 x 10~ n 7 



B 



(25) 



(26) 



,77 AXJJ ' V 10 ~ 1G , 

which gives, roughly, a gyroradius (in units of the plasmon 
size) of 10 -9 . The diffusive escape time with such a mean 
free path is then 

R 2 



to 



D( 7 ) 
(lO 2 ") 



R 



77 AU 



B 



lO-iG 



(27) 



10 years. 

Spectral steepening cannot be due to diffusive escape 
if the scattering occurs on the scale of a gyroradius. A 
leakage timescale of a few days would require a mean free 
path of 

A - IQ^R, (28) 

thus leakage via this mechanism is not a plausible expla- 
nation for the observed energy losses, given the gyroradius 
derived in Equation 26. 

9.5. Light travel-time effects 

We now consider an alternative explanation for the ob- 
served high-frequency spectral evolution, namely that it 
is governed by light travel-time effects in a series of plas- 
mons which are ejected and radiate at different distances 
from the observer, and evolve from being optically-thick 
to optically-thin as they move outwards from the core and 
expand (in fact the effect will be important whether or 
not the plasmons expand, as long as they have a finite 
size and get fainter with time. We invoke expansion, since 
that is implied by our observations, as discussed in § 6.3.) 
At 5 GHz (the frequency at which we have best temporal 
coverage of VLBA images) , the separation of the northern 
and southern emitting regions expands from a projected 
size of 12 mas on September 18 to 39 mas on September 
21. For a source distance of lOkpc, this equates to 120 AU 
on September 18 and 390 AU on September 21. The sizes 
of the individual plasmons increase from 7 mas to 17 mas, 
equating to 70 and 170 AU respectively. Thus it takes 
light a minimum of between 9 and 22 hours respectively to 
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travel across each plasmon, and between 16 and 52 hours 
to travel between the plasmons. The source is evolving 
on a similar timescale, so light travel-time effects must be 
relevant. Since the inclination of the jet axis to the line of 
sight has been found to be small (i = 10.5 ± 4.2°, § 6.4.3), 
light travel-time effects between the plasmons will be im- 
portant. At any one time, we see light from the far side of a 
given plasmon at an earlier epoch than light emitted from 
its near side. At earlier epochs the plasmon's magnetic 
field and emissivity were higher. As the plasmon expands, 
the field, emissivity, and any opacity decrease. If the com- 
bination of size and magnetic field are such that opacity 
effects are initially important, then initially an overall flat 
(or inverted) spectrum would be seen. During the tran- 
sition from t >• 1 to t < 1, the near side of a plasmon 
will contribute optically thin emission, but we will only 
see in to an optical depth of r « 1. After sufficient time 
has elapsed, when all of the plasmons are optically thin, 
a constant, steep, terminal value of the spectral index will 
be observed. 

Our toy model is shown in Figure 9; our analysis follows 
that of Blundell & Alexander (1994) who considered anal- 
ogous effects in the evolution of the lobes of radio galaxies. 
We consider light setting off from the closest point of the 
first plasmon to the observer, and determine what time 
the light would have left other parts of that plasmon, and 
other plasmons, in order to reach plane / at the same time. 
All light leaving plane / (perpendicular to the line of sight) 
at a given instant clearly reaches the observer simultane- 
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Fig. 9. — Assumed geometry of the jet. C is the core of the 
system, and the observer is in direction O. The plasmon centres 
move outwards at speed f)c and the jet axis is aligned at angle i to 
the line of sight. We use the co-ordinate y as a dimensionless local 
measure of the distance from the front surface to any plane within a 
given plasmon. Light from position y within a plasmon leaves that 
plasmon in the approaching jet at time t ya in order to reach plane / 
at time tf . Light from position y within the corresponding plasmon 
in the receding jet leaves at an earlier time t yT in order to reach / 
at this same time. All light leaving / at time tf will obviously be 
received at the same time by the observer. 



ously. The refractive index of the plasmons is taken to be 1 
(in reality, it will be greater than 1, and this will make the 
light travel-time effects even more significant). We assume 
a tangled field, so that the magnetic flux density falls off 
as B cx R^ 1 (as described in §7.2.1). Given the magnetic 
field at any one time, and having calculated the time the 
light set off from a plane within a plasmon, and knowing 
the time at which the plasmon was ejected, then assuming 
constant uniform expansion of the plasmons, the size of a 
plasmon at any given time, and therefore its magnetic field 
can be determined. Dividing each plasmon up into slices 
perpendicular to the line of sight, from which light sets off 
progressively earlier, we derive the field in each slice at the 
time at which light leaves that slice. We can then calcu- 
late the emissivity of each slice, assuming the plasmons to 
be homogeneously filled with energetic particles. We then 
integrate over the entire plasmon to find the flux density 
expected at any particular frequency, having corrected for 
Doppler boosting and relativistic time-dilation effects. 

To model the emissivity of each slice, rather than take 
the full expression for the synchrotron emissivity, we as- 
sume that each electron radiates at a single frequency, 
determined by its Lorentz factor and the ambient mag- 
netic field (Equation 9). In order to model the opacity, 
we put in an absorption coefficient (derived either for syn- 
chrotron self-absorption, or for free-free absorption, as in 
Rybicki & Lightman (1979)). We calculate the location of 
the r = 1 surface within each plasmon at each frequency, 
and assume we see no emission from beyond that surface. 
Emission from the near side of the r = 1 surface is taken 
to be synchrotron emission from a power-law distribution 
of electrons with electron index p = 2.2. This produces 
a turned-over spectrum as we observe. We then observe 
how the spectrum evolves with time. 

Plots of spectra generated are given in Figure 10. The 
full theory will be explored in more detail in a companion 
paper (Miller- Jones et al. in prep.). For now, we note that 
such a model produces the qualitative features of spectra 
similar to those observed, in that it produces a turnover 
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Fig. 10.— Model spectra for 2001 September 18, 19, 21, 22, and 
23, to be compared with Figure 6. Spectral indices between each 
sampled frequency point are given. The model takes into account 
light travel-time and opacity effects for a series of jet knots, with 
adiabatic expansion of a scries of ejected plasmons governing the 
spectral evolution. As shown above, it does not include any flat- 
spectrum core component. 



Radio Observations of Cygnus X-3 in Flare 



19 



moving to lower frequencies with time and a steepening 
of the spectrum at frequencies just above the turnover to 
a terminal value of the spectral index which corresponds 
to the underlying electron energy index. An important 
success of this model is that the timescales for the move- 
ment of the turnover and for the steepening of the high- 
frequency spectrum match those observed. 

However, where this model falls short is that it does 
not predict a spectrum flatter than v~ 6 all the way up 
to the highest observing frequencies. However, this model 
only includes jet knots which evolve as described above: it 
takes no account of any flat-spectrum core component. We 
note that invoking the addition of a core with a canonical 
v~ Q1 spectrum whose intensity decreases with time might 
sufficiently flatten the predicted high-frequency spectrum 
at early epochs. Such a core component was postulated 
in §6.4.4 (see also Figure 5) and was present in this out- 
burst, as we now describe. In §6.5, we found that the ob- 
served VLBA knots appeared to have been ejected 2.5±0.5 
days after the integrated radio flux density began to rise. 
This would be consistent with a scenario in which a flat- 
spectrum core component as postulated above dominated 
the flux density initially, but faded over time. After the 
ejection of the steep-spectrum radio plasmons, their flux 
density would later come to dominate the observed flux 
density, especially at lower frequencies. 

10. COMPARISON TO RADIO GALAXIES 

We now make a brief comparison between the magnetic 
field and energy budget in this microquasar flare and those 
parameters in quasars: scaled-up, longer-lived versions of 
Cygnus X-3. 

10.1. Magnetic Field 

M 87 is a radio galaxy at a distance of 16Mpc for which 
there arc good constraints on the magnetic field at different 
distances from the core, which has a central black hole of 
mass 3.2 ±0.9 x 10 9 M Q (Macchetto et al. 1997). Reynolds 
et al. (1996) derived limits to the magnetic field in the jet 
of M87 of 0.01 < B < 0.2 G at an angular separation of 
0.1 mas from the core, approximately 1800 AU, or 40 gravi- 
tational radii (re) from the central engine. Heinz & Begel- 
man (1997) derived a magnetic field of 5 < B < 49 y(iG 
at several kpc (<~ 10 6 — 10 7 rc) from the core. Standard 
cquipartition arguments give 300 < B < 600 ^G at that 
distance (e.g. Meisenheimer et al. 1996). 

Cygnus X-3 is thought to contain either a neutron star or 
a black hole of a few solar masses (e.g. van den Heuvel & de 
Loore 1973; Schmutz et al. 1996; Ergma & Yungelson 1998; 
Stark & Saia 2003): a factor of 10 8 - 10 9 times smaller 
than the black hole in M 87. Assuming that the compact 
object in Cygnus X-3 has a mass of a few solar masses, 
its gravitational radius is of order 10 km. Observing with 
a few mas resolution probes the source at 10 8 — 10 9 rQ. 
According to Heinz (2002), magnetic fields in the inner jet 
and disk should scale as M -1 / 2 where M is the black hole 
mass. Comparing the magnetic fields in the inner jets of 
M87 and Cygnus X-3 this predicts the field in the inner 
jet (tens of gravitational radii from the compact object) 
of Cygnus X-3 to be of order 10 3 G. This is a factor of 10 4 
higher than that which we observe further out at <~ 10 8 tq- 
Let R be the jet cross-section and r the distance from the 



core measured along the jet axis. Taking a jet model in 
which the cross-section scales as R cx r 5 , and invoking flux- 
freezing and a tangled magnetic field, such that B oc R~ x , 
we find that B oc r~ s . Putting in the values we have for 
B and r, we find that 0.50 < S < 0.83. This limit may 
be compared with theoretical limits quoted by Heinz & 
Sunyaev (2002). They found that 

R = R (r/r f, < f3 < 1 (29) 

n = n (r/r )- x , 2/3 < A < 3/3 (30) 

B = B a (r/r a y S , 1/2A < 8 < A (31) 

where they defined r as the distance from the black hole, r 
as the distance at which the jet is injected (taken as 10 tq) , 
R as the jet cross-section and n as the particle number 
density. They did however point out that this assumes self- 
similarity. It would thus be expected to hold in the inner 
jet, but not necessarily further out where the jet behaviour 
may be dominated by interactions with its environment. 
The case [3 = 1 corresponds to a jet with constant opening 
angle. A jet composed of discrete ejections has A ~ 3/3, 
whereas a continuous jet would satisfy A ~ 2(3. A value 5 = 
2/3 A is predicted for a completely tangled magnetic field. 
From the VLBA observations in Section 5, we observed a 
jet which seems to be composed of discrete ejections, and 
we have assumed a tangled field. Thus we would expect 
S = 2/3. Our derived limits on the value of the index S 
then imply 0.25 < (3 < 0.41. Since the VLBA images did 
not have sufficient transverse resolution to monitor the jet 
opening angle, and since they showed a series of discrete 
knots in the jet, it is not possible to verify accurately any 
variation in the opening angle from the observations. 

We note, however, that although we expect similar physics 
to apply in the inner, jet-producing regions of AGN and 
microquasars, the environments of the two types of sys- 
tem arc different. According to Heinz (2002), in dynam- 
ical terms, microquasar jets are located in relatively low- 
pressure, low-density environments. Where we observe the 
radio jets (i.e. as they interact with their environments), 
the characteristics of the environment will affect the kine- 
matics of the jet, and a simple scaling between the two 
types of system is not expected. 

Heinz & Begelman (1997) found that the magnetic field 
in the M 87 jet at a distance of several kpc from the core 
should be a factor ~ 1.5 — 5 times lower than the equipar- 
tition value. We found the Cygnus X-3 outer jet to be out 
of equipartition (modulo the uncertainties in this calcula- 
tion) by a similar factor of a few (§8.4). 

10.2. Energetics 

The typical energy thought to be stored in the lobes 
of a radio galaxy such as Cygnus A is estimated to be of 
order 10 59 erg (Kembhavi & Narlikar 1999, chapter 9), and 
the typical kinetic jet power of order 10 42 — lO^ergs" 1 
(Reynolds et al. 1996). From the results of September 18, 
it is possible to deduce the minimum energy present in 
the outburst, using the 7 max and 7 m ; n corresponding to 
the highest and lowest observed frequencies on that date. 
This gives a lower limit of 8.5 x 10 42 erg for the energy 
released following this flare. We do not believe, however, 
that we have fully sampled the range of Lorentz factors 
present in the jet, so this should be taken as a lower limit, 
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and fuller sampling would be required in order to make a 
complete comparison with AGN. 

11. CONCLUSIONS 

We have presented time-sequenced milliarcsecond-resolution 
observations at three different frequencies of a major out- 
burst of Cygnus X-3, accompanied by simultaneous quasi- 
daily monitoring of the (unresolved) flux density of the 
source, with unprecedented spectral coverage of a large 
flare, ranging over three decades in frequency. The proper 
motions of individual plasmons have been directly tracked 
for the first time, and the expansion of those plasmons has 
been directly measured. 

One main result is to securely determine the presence 
of a two-sided jet in Cygnus X-3 on milliarcsecond scales 
following a major flare. The southern jet appears to be ap- 
proaching and the northern counterjet appears to be reced- 
ing from us, being substantially de-boosted. We derived 
a jet speed of (3 ~ 0.63. The jet appears to be precessing 
with a period of order 5.3 days, significantly shorter than 
that derived by Mioduszewski et al. (2001) for the 1997 
flare. It appears that the milliarcsecond-scale jet is not 
ejected immediately on the rise of the integrated flux den- 
sity, but ^2.5 days after this, a point which we explore in 
a forthcoming paper (Miller- Jones et al. in prep.). 

We believe that it is possible that the intrinsic character- 
istics of Cygnus X-3 in the outbursts of, for example, 1997 
and 2001, may have been different. It is not clear that all 
outbursts of such sources follow the same patterns. Cer- 
tainly the flux densities and fading times are not the same 
for different events, as outlined above. This variable, al- 
most whimsical, nature of the source needs to be clarified. 
Certain parameters would not be expected to change (e.g. 
after accounting for precession, the inclination angle of the 
jet axis to the line of sight might be expected to be con- 
stant). Identical spectral, temporal, and «u-coverage of 
different events with similar radio lightcurves might help 
to determine whether the jet parameters are indeed con- 
stant for different flares. 

Observations of the integrated flux density of the source 
showed an evolving spectrum similar to that seen in pre- 
vious outbursts of Cygnus X-3. The high-frequency spec- 
trum at frequency v took the form v~ a , where a steepened 
from <~ 0.4 to ~ 0.6 and the spectral index then remained 
at this latter terminal value. We contend that the spectral 
evolution is the result of light travel-time effects in plas- 
mons which were initially optically thick, but which be- 
came optically thin as they expanded. Such a model can 
produce spectral evolution over timescales similar to those 
observed. Confirmation that opacity effects are important 
came from the integrated spectrum between 1.4 GHz and 
74 MHz: this showed a turnover moving to lower frequen- 
cies with time, caused by the decreasing opacity (probably 
caused by either synchrotron self-absorption, free-free ab- 
sorption or a low-energy cutoff in the electron spectrum). 
The presence of an additional flat-spectrum core compo- 
nent which was responsible for the initial rise of the flare 
but which then faded with time was also postulated. 

We derived an upper limit to the magnetic field in the 
plasmons of 0.15 ± 0.06 G, consistent with previous esti- 
mates for similar systems. We also placed lower limits on 
the minimum and maximum observed Lorentz factors of 



the emitting electrons, of 7 = 58 and 7 = 321. 

Despite the unique nature of this system, Cygnus X-3 
does share many of the characteristics of other, better un- 
derstood microquasars. Owing to its position within the 
Galaxy and the associated difficulties in imaging the sys- 
tem, Cygnus X-3 still remains somewhat of an enigma. 
Studying it can, however, enhance our understanding of 
the effects which we believe must play a part in the ob- 
served emission from microquasars. Light travel-time ef- 
fects are unavoidable in systems of this size evolving on 
similar timescales. 
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Table 1 



Cygnus X-3 VLA Flux Densities 



Date (UT) 




JD-2450000 


Observing Frequency (GHz) 


Flux Density (Jy) 


2001 September 


09 


2161.73955 ± 0.00050 


4.860 


0.0209 ± 0.0007 






21ol. ( 0000 ± 0.0004b 


8.460 


O.Oolz ± O.OOOo 






2161.73212 ± 0.00075 


14.94 


r\ ft A ft a t ft nfttft 

0.0424 ± 0.0020 






2151.72419 ± 0.00127 


22.46 


n 1 ft nfftft 

0.0387 ± 0.0030 






2161.71748 ± 0.00139 


43.34 


0.0472 ± 0.0038 


2001 September 


18 


2170.57112 ± 0.00168 


1.385 


8.56 ± 0.17 






2170.57112 ± 0.00168 


1.465 


9.16 ± 0.18 






2170.57570 ± 0.00093 


4.860 


13.41 ± 0.27 






01 7n c^777q 4- n nnnQ 1 ? 


0.40U 


11.00 ± u. Iz 






2170.59215 ± 0.00091 


14.94 


9.04 ± 0.43 






2170.58810 ± 0.00197 


22.46 


7.44 ± 0.57 b 






2170.58351 ± 0.00168 


43.34 


5.99 ± 0.46 


2001 September 


19 


2171.45749 ± 0.00050 


1.385 


14.25 ± 0.28 






2171.45749 ± 0.00050 


1.465 


14.31 ± 0.29 






2171.45660 ± 0.00058 


4.860 


13.43 ± 0.27 






2171.46615 ± 0.00052 


8.460 


10.54 ± 0.11 






2171.48380 ± 0.00058 


14.94 


8.22 ± 0.39 






2171.48079 ± 0.00127 


22.46 


6.70 ± 0.51 b 






2171.47454 ± 0.00127 


43.34 


4.62 ± 0.37 


2001 September 


21 


2173.63970 ± 0.00058 


1.385 


7.98 ± 0.16 






2173.63970 ± 0.00058 


1.465 


7.86 ± 0.16 






2173.64092 ± 0.00041 


4.860 


4.83 ± 0.10 






2173.65000 ± 0.00046 


8.460 


3.54 ± 0.04 






2173.66476 ± 0.00041 


14.94 


2.37± 0.11 






2173.66119 ± 0.00050 


22.46 


1.79 ± 0.14 






2173.65567 ± 0.00046 


43.34 


1.22 ± 0.15 


2001 September 


22 


2174.52118 ±0.01806 


0.3275 


6.161 ± 1.4 






2174.50232 ± 0.00058 


1.385 


5.640 ± 0.113 






2174.50232 ± 0.00058 


1.465 


5.377 ± 0.108 






2174.52558 ± 0.01782 


4.860 


2.968 ± 0.059 






2174.53115 ± 0.01587 


8.460 


2.126 ± 0.021 






2174.53671 ± 0.00091 


14.94 


1.465 ± 0.070 






2174.52646 ± 0.00108 


43.34 


0.688 ± 0.055 


onm — r 1 — i r — 

2001 September 


23 


2175.70995 ± 0.01053 


n n7QO 

0.0738 


< 2.7" 






2175.70723 ± 0.01048 


0.3275 


4.061 ± 1.4 






2175.69601 ± 0.00052 


1.385 


2.919 ± 0.058 






2175.69601 ± 0.00052 


1.465 


2.834 ± 0.057 






2175.70139 ± 0.00046 


4.860 


1.425 ± 0.029 






2175.70642 ± 0.00041 


8.460 


0.990 ± 0.010 






2175.71570 ± 0.00062 


14.94 


0.641 ± 0.031 






2175.71250 ± 0.00058 


22.46 


0.492 ± 0.038 






2175.71078 ± 0.00091 


43.34 


0.315 ± 0.030 


2001 September 


24 


2176.69213 ± 0.01563 


0.0738 


< 4.5" 






2176.68918 ± 0.01557 


0.3275 


7.001 ± 1.4 






2176.69371 ± 0.02157 


1.385 


1.893 ± 0.038 






2176.69371 ± 0.02157 


1.465 


1.842 ± 0.037 






2176.69618 ± 0.01794 


4.860 


0.895 ± 0.018 






2176.69699 ± 0.03117 


8.460 


0.638 ± 0.006 






2176.69884 ± 0.00058 


14.94 


0.470 ± 0.022 






2176.69220 ± 0.00065 


22.46 


0.392 ± 0.030 






2176.69041 ± 0.00091 


43.34 


0.314 ± 0.016 



a Thcsc are 3(7 upper limits as no source was detected at the known position of Cygnus X-3 
b Tipping scans were used to derive the atmospheric opacities for these observations 



Note. — The error bars given were derived as described in 3.3 



Table 2 

Cygnus X-3 OVRO Millimeter. Array Flux Densities at 99.48 GHz 



Date (UT) JD-2450000 Flux Density (Jy) 

2001 Sep 20 2172.6930 ± 0.0597 1.206 ± 0.005 

2001 Sep 21 2173.6242 ± 0.0167 0.942 ± 0.007 

2001 Sep 22 2174.8232 ± 0.0170 0.363 ± 0.010 

2001 Sep 23 2175.6892 ± 0.0205 0.216 ± 0.005 

2001 Sep 26 2178.5559 ± 0.0275 0.121 ±0.004 



Note. — The uncertainties quoted in the flux density are 
those due to thermal noise. 



Radio Observations of Cygnus X-3 in Flare 



23 



Table 3 



VLBA IMAGE PARAMETERS 



Date (UT) 


JD-2450000 


v (GHz) 


Beam (mas x mas) 


Beam P.A. 


Scattering Disk Size (mas) 


r.m.s. (mjy bm 1 ) 


2001 Sop 18 


2170.78779 ± 0.17781 


22.233 


8.59 x 3.48 


-59.4 


0.69 ± 0.22 


19.7 


2170.76849 ± 0.18423 


4.995 


13.71 x 12.92 


+ 16.7 


15.54 ± 2.42 


13.8 


2001 Sep 19 


2171.61065 ± 0.20914 


22.233 


8.30 x 4.00 


-29.1 


0.69 ± 0.22 


9.68 


2171.59074 ± 0.20183 


15.365 


16.19 X 8.76 


-19.8 


1.48 ± 0.42 


126.4 




2171.59149 ± 0.21542 


4.995 


15.65 x 13.58 


-63.9 


15.54 ± 2.42 


6.38 


2001 Sop 20 


2172.75584 ±0.20548 


22.233 


5.88 x 3.77 


-59.2 


0.69 ±0.22 


6.52 




2172.72633 ± 0.18880 


15.365 


9.39 x 8.91 


+ 19.2 


1.48 ± 0.42 


42.0 




2172.71344 ± 0.18880 


4.995 


14.54 x 13.45 


-85.8 


15.54 ± 2.42 


17.9 


2001 Sop 21 


2173.65190 ± 0.26304 


15.365 


13.31 x 9.00 


-14.4 


1.48 ±0.42 


21.0 




2173.63907 ± 0.26305 


4.995 


11.85 x 10.60 


+61.4 


15.54 ± 2.42 


7.49 



Table 4 



Fitted Gaussian FWHM sizes of the 22 GHz VLBA knots 



Date 


Component 


Fitted FWHM (mas)' 


Major Axis (mas) 


Minor Axis (mas) 


P.A. 


2001 Sep 18 
2001 Sep 19 
2001 Sep 19 
2001 Sep 20 
2001 Sep 20 


Single component 
North 
South 
North 
South 


7.88 ± 0.51 
11.98 ± 0.15 

13.36 ± 0.21 

14.37 ± 0.54 
16.92 ± 0.88 


7.23 ± 0.17 
11.00 ±0.14 
10.93 ± 0.23 
10.58 ± 0.39 
13.89 ± 0.71 


2.65 ± 0.34 

4.03 ± 0.09 

5.04 ± 0.12 
6.02 ± 0.39 
6.36 ± 0.58 


167.5 ± 2.0 
168.1 ± 0.7 
157.5 ± 1.2 
26.5 ± 94.3 
25.0 ± 95.8 



a The diameter of a sphere subtending the same solid angle as the Gaussian with the given major and minor axes, 
calculated according to Equation 8. 



Note. — Deconvolved sizes and errors were derived using the AIPS task JMFIT 



Table 5 

Results of fitting the 22 GHz VLBA image of September 20 with the Hjellming & Johnston processing jet model 



Core Ft. A. 


Core Dec. 


Srot 


Parameter 


Mean 


Standard Deviation 


Maximum 


Minimum 


20 32 25.77335 


+40 57 27.9650 


-1 




168.8 


91.1 


345.9 


3.3 








■>!> 


2.4 


1.2 


5.8 


0.4 








i 


10.5 


4.2 


18.8 


2.5 








X 


98.9 


2.9 


106.2 


93.4 








13 


0.63 


0.09 


0.80 


0.47 








P (days) 


5.34 


0.67 


7.15 


3.73 








d (kpc) 


10.8 


2.6 


14.8 


5.6 








tview (days) 


7.44 


0.54 


8.43 


6.08 


20 32 25.77335 


+40 57 27.9650 


+ 1 


4> 


136.1 


97.7 


330.2 


2.1 








V- 


2.2 


0.9 


4.3 


0.8 








i 


169.6 


3.8 


175.6 


161.3 








X 


278.1 


2.8 


283.6 


273.6 








& 


0.60 


0.08 


0.77 


0.49 








P (days) 


5.49 


0.46 


7.06 


4.81 








d (kpc) 


10.0 


2.1 


13.0 


6.2 








t viow (days) 


7.53 


0.52 


8.42 


6.61 



Note. — All angles (0, ip, i, and x) are measured in degrees 



